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I.  SUMMARY 

This  report  describes  the  operation,  maintenance  ami  research  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  October  1984  -  31 
March  1985. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged 
98.2  per  cent  as  compared  to  95.1  for  the  previous  period.  Disk  failure  and 
maintenance  in  connection  with  disk  failure  caused  most  of  the  downtime. 

A  total  of  1792  events  were  reported  in  this  period,  giving  a 
daily  average  of  9.9  events. 

During  this  period  the  4331  IBM  computer  has  been  replaced  by  a  more 
powerful  4341.  The  increased  capacity  will  enable  us  to  improve  the  pro¬ 
cessing  algorithms,  test  processing  together  with  online  operation  and  will 
make  program  development  easier.  The  new  automatic  bulletin  transfer  has 
been  running  for  approximately  8  months,  and  experience  so  far  has  been 
good. 

The  majority  of  maintenance  visits  to  the  subarrays  have  been  in  connection 
with  corrective  maintenance,  including  repair  of  cables. 

The  research  activity  is  summarized  in  section  VII.  Section  VII. 1 
discusses  real  time  event  detection  using  the  small-aperture  NORESS  array, 
while  VII. 2  presents  a  preliminary  evaluation  of  the  event  detection  and 
location  capabilities  of  NORESS.  Section  VII. 3  takes  up  seismic  modelling 
of  an  anisotropic  boundary  layer,  and  Section  VII .4  discusses  an  application 
of  long-period  waveform  inversion  to  shallow  sources  under  the  Tibetan 
Plateau.  In  section  VII. 5  an  Interactive  analysis  program  for  use  with 
NORESS  data  is  presented.  Analysis  of  teleseismic  P  coda  using  NORSAR  and 
NORESS  is  discussed  in  section  VII. 6.  Results  from  3-component  seismogram 
analysis  are  presented  in  section  VII. 7,  while  epicenter  location  scenarios 
are  presented  in  section  VII. 8.  In  section  VII. 9  the  current  status  of  the 
RST/RSTE  research  project  is  presented.  A  study  of  the  crustal  structure  in 
the  NORSAR  array  siting  area  derived  from  gravity  observations  is  described 
in  section  VII. 10. 


-  2  - 


II.  OPERATION  OF  ALL  SYSTEMS 

II.  1 _ Detection  Processor  (DP)  Operation 

There  have  been  125  breaks  in  the  otherwise  continuous  operation  of 
the  NORSAR  online  system  within  the  current  6-month  reporting 
interval.  Disk  failure  and  maintenance  in  connection  with  disk  failure 
caused  most  of  the  downtime.  The  uptime  percentage  for  the  period  is 
98.2  as  compared  to  95.1  for  the  previous  period. 

Fig.  II. 1.1  and  the  accompanying  Table  II. 1.1  both  show  the  daily  DP 
downtime  for  the  days  between  1  October  1984  and  31  March  1985.  The 
monthly  recording  times  and  percentages  are  given  in  Table  11.1.2. 

The  breaks  can  be  grouped  as  follows: 


a)  MODCOMP  failure  :  4 

b)  Disk  failure  :  31 

c)  Maintenance  stops  :  16 

d)  Power  jumps  and  breaks  :  3 

e)  Stops  related  to 

system  operation  :  4 

f)  TOD  correction  :  49 

g)  Communications  lines  :  18 


The  total  downtime  for  the  period  was  77  hours  and  30  minutes.  The 
mean-time-between-f ailure  (MTBF)  was  1.4,  the  same  as  for  the  previous 
period. 
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Month 

DP  Uptime 
(hrs) 

DP  Uptime 
<%) 

No.  of 

DP  breaks 

No.  of  days 
with  breaks 

DP 

MTBF* 

(days) 

Oct  84 

740.05 

99.5 

20 

14 

1.5 

Nov 

717.12 

99.6 

15 

11 

1.9 

Dec 

743.44 

99.9 

11 

11 

2.6 

Jan  85 

717.54 

96.5 

34 

18 

0.9 

Feb 

652.01 

97.0 

20 

11 

1.3 

Mar 

719.34 

96.7 

25 

12 

1.2 

4290.30 

98.2 

125 

77 

1.4 

Mean-time-between-failures  =  (total  uptlme/no.  of  up  intervals) 


Table  II. 1.2  Online  system  performance,  1  October  1984  -  31  March 


-  7  - 


II. 2 _ Array  Communications 

The  performance  of  Che  communications  systems  throughout  the  period  is 
reflected  in  Table  11.2.1.  As  the  table  indicates,  the  majority  of  the 
systems  have  been  affected,  mainly  by  broken  cables  or  defective  com- 
pon<  _s  causing  unwanted  effects  such  as  delays  and  attenuation.  Other 
incidents  have  been  caused  by  NTA,  mainly  in  connection  with  replace¬ 
ment/reconfiguration  of  equipment  situated  in  the  telestations  at 
Lillestr0m,  Hamar  and  Lillehammer.  Tests  and  changes  at  the  NDPC  have 
also  affected  the  communications  system's  uptime.  Systems  marked  with 
asterisks  are  commented  in  the  following  paragraphs. 

Summary 

Oct.  84  Most  systems  were  affected,  02B,  02C  and  04C  significantly. 

Also  01A,  01B  and  03C  had  outages.  Degradation  of  the  02B  and 
04C  line  performance  started  in  August/September  84,  and  in 
both  cases  damaged  pupin  coils  caused  mismatch  problems. 

Tests  indicated  improved  conditions  for  both  systems  in  C- 
loop  tests  (digital  loop)  15  and  23  oct.,  respectively. 
Therefore  the  systems  were  temporarily  reconnected  awaiting 
NTA/ Hamar' s  next  initiative.  02C  was  disconnected  from  the 
MODCOMP  30  October  due  to  excessive  error  rate  caused  by  work 
going  on  at  the  NTA  Lillehammer.  Otherwise,  scheduled  work  at 
NTA/Lillestrpm  affected  all  communications  systems  weeks  41 
and  42.  06C  ordinary  comm,  line  section  between  Hamar  and  the 
NTA  central  near  Rokosjpen  was  by  a  mistake  disconnected  in 
the  beginning  of  August  1984,  when  not  needed  for  transfer  of 
analog  data  from  the  NORESS  experimental  array.  When  starting 
recording  from  06C  (NORSAR)  again,  the  short  line  between  the 
CTV  and  Rokosj0en  central  was  connected  to  the  PCM  equipment, 
which  in  the  meantime  had  been  installed  there.  Afterwards 


NTA  has  not  been  able  to  establish  an  acceptable  com¬ 
munications  channel.  03C  error  figures  were  caused  by  discon¬ 
nection  of  the  M0DC0MP  from  the  communications  system  for 
test  purposes. 

Nov.  84  Also  this  period  NTA/Lillehammer  caused  02C  problems.  Spikes 
were  observed  5  November  and  communications  tests  proved  bad 
performance  in  C-loop  (digital  loop  on  CTV  modem).  6  November 
point-to-point  test  (Kjeller-Hornsj0)  was  carried  out.  Modems 
at  both  locations  found  in  order;  the  line  towards  the  CTV 
was  suspected.  In  connection  with  work  at  Lillehammer 
telestation,  including  change  of  channel  equipment,  the  ori¬ 
ginal  equalizers  had  not  been  reinstalled.  The  subarray  was 
back  in  operation  7  November. 

02B.  Spikes  observed  7  November.  NTA/Hamar  admitted  last  sec¬ 
tion  of  cable  remained  to  be  check.  Fault  location  started 
but  not  finished  (8  November).  However,  test  proved  notice¬ 
able  improvement,  but  still  a  few  errors  in  C-loop.  The 
subarray  was  demasked  and  back  in  operation.  NTA/Hamar 
resumed  work  12  November  and  a  "blown"  pupin  coil  was  found 
on  the  cable  section  near  the  CTV.  Still  a  few  errors  were 
observed  in  the  C-loop. 

01B.  13  November  spikes  were  observed  on  all  channels  due  to 
a  damaged  cable,  also  input  level  to  Kjeller  low  (-  35  dB). 

02B.  Communications  cable  broken  15  November.  Cable  repaired 
and  subarray  back  in  operation  the  same  day.  However,  spikes 
in  data  were  again  observed  on  all  channels  (SP).  NTA/Hamar 
found  3  more  "blown"  pupin  coils  on  the  last  cable  section. 
Finally,  22  November  test  proved  error-free  system. 

02C.  The  system  was  down  for  about  2.5  hours  30  November.  The 
same  day  all  systems  were  down  for  a  while.  Afterwa  Is  04C 
did  not  resume  operation.  Modem-to-modem  test  stated  error- 
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free  communications  system.  The  conclusion  was  that  a  SLEM 
(subarray  electronic  module)  "hang-up"  had  occurred.  After  a 
restart  of  the  communications  machine  (MODCOMP),  the  system 
was  operational  again. 


Dec.  84  01B  resumed  operation  20  December  after  having  been  down 

since  14  November.  All  communications  systems  were  affected 
28  December  due  to  CF-channel  change  at  Lillestr0m  and  Hamar 
telestations.  06C  was  still  out  of  operation. 


Jan.  85  01A.  In  the  period  4-8  Jan  we  were  not  able  to  establish  con¬ 

tact  with  the  subarray.  NMC  staff  visited  the  subarray  8 
January  and  found  in-level  reduced  by  nearly  50%.  A  faulty 
equalizer  was  found  and  replaced  (NTA/Hamar). 

01B.  15  January  the  level  toward  the  subarray  was  increased 
(on  CF-channel  level)  as  spikes  in  data  were  observed. 

06C.  As  an  attempt  to  solve  the  06C  communications  problem, 
new  NOKIA  2400  baud  modems  were  installed  9  January  (NDPC  and 
CTV) .  After  the  installation  a  test  was  run  modem-to-modem, 
apparently  with  acceptable  results.  When  reconnecting  the 
NDPC  and  CTV  data  equipment,  error  figures  far  outside  speci¬ 
fications  were  observed. 


Feb.  85  06C.  In  order  to  eliminate  as  many  error  sources  as  possible, 

the  06C  SLEM  was  replaced  27  February.  As  no  change  in  the 
error  pattern  was  observed,  the  original  unit  was 
reinstalled. 


Otherwise  high  performance  of  all  systems  In  February. 

Mar.  85  02C.  Scheduled  system  work  at  Llllehammer  telestation  caused 

unwanted  outages  several  times  this  period.  Between  11  and  15 
March  the  02C  communications  system  remained  down.  19  and  28 
March  stops  of  approx.  2  hours  duration  were  observed.  A 
shorter  outage  occurred  22  March. 

01A  communication  which  stopped  19  March  was  back  in  opera¬ 
tion  28  March.  A  bad  cable  caused  this  outage. 

02B.  Line  work  in  the  Hamar  area  caused  an  outage  21  March. 
02B  (telemetry)  stations  resumed  operation  22  March  after 
having  been  down  since  18  February  due  to  a  bad  com¬ 
munications  cable.  The  reason  for  why  the  cable  repair  was 
postponed  had  connection  with  access  problems  caused  by  deep 
snow  in  the  area. 

06C.  NTA/Lillestr0m  and  Hamar  strapped  the  two  new  NOKIA  2400 
modems  to  external  clock,  which  is  correct.  We  therefore 
believed  that  problem  now  was  solved.  The  performance  of  the 
communications  system  is  even  worse  than  before  the 
strapping.  We  have  therefore  in  a  letter  to  NTA  recently 
asked  them  to  reestablish  the  original  communications  path 
between  Hamar  and  Rokoberget  central. 

O.A.  Hansen 


III.  ARRAY  PERFORMANCE 

III.l  Event  Processor  Operation 

In  Table  III.l  some  monthly  statistics  of  the  Event  Processor  opera 
tion  are  given: 


Teleseismic 

Core  Phases 

Sum 

Daily 

Oct  84 

216 

54 

270 

8.7 

Nov  84 

196 

66 

262 

8.7 

Dec  84 

356 

53 

409 

13.2 

Jan  85 

214 

44 

258 

8.3 

Feb  85 

215 

52 

267 

9.5 

Mar  85 

271 

55 

326 

10.5 

1468 

324 

1792 

9.9 

Table 

111.1 

Event 

Processor  Statistics,  Oct  1984 

-  Mar 

1985 

B.  Paulsen 


IV.  IMPROVEMENTS  AND  MODIFICATIONS 

IV.  1 _ NORSAR  Online  System  using  4331/4341  and  MODCOMP  Classic 

We  refer  to  the  detection  processor  operation  statistics  for  detailed 
information  about  uptime,  communication  and  processing  performance. 

Up  to  January  1985  the  NORSAR  online  and  detection  system  consisted  of 
oe  IBM  4331  processing  computer  with  a  high-speed  link  to  a  front  end 
communication  processor  which  received  data  from  all  SLEMs.  During 
this  report  period  the  4331  has  been  replaced  by  a  more  powerful  com¬ 
puter,  the  IBM  4341.  The  increased  capacity  gives  us  the  opportunity 
to  improve  the  processing  algorithms,  test  processing  together  with 
online  operation  and  makes  program  development  easier.  Together  with 
the  new  processing  system,  we  have  the  following  main  new  features  for 
the  NORSAR  online  system: 

-  new  disc  controller 

-  data  from  06C  reestablished  (06C  previously  used  for  NORESS 
developments) 

improved  alarm  system 

-  experience  with  new  bulletin  transfer 
ARPANET  connection. 

New  computer 

The  4331  has  been  replaced  with  our  previous  research  computer  (4341). 
The  upgrade  gives  us  the  opportunity  to  more  easily  change  the  NORSAR 
online  system  with  a  minimum  of  processing  stops.  The  4341  is  about  4 
times  the  4331  computer  in  CPU  speed.  The  CPU  speed  is  close  to 
1  Mips.  The  memory  size  on  the  IBM  4341  computer  is  4  Mb  while  the 
4331  only  had  1  Mb.  The  increase  memory  size  makes  it  easier  to  run 
additional  programs  together  with  the  online  processing. 


New  disc  controller 

In  the  4331  system  the  disc  controller  was  built  into  the  computer. 

The  built-in  disc  controller  had  some  limitations  such  as  addressing 
range,  speed  and  high  load  on  CPU.  Together  with  the  new  computer  (IBM 
4341)  we  installed  a  new  separate  disc  controller  (3694)  from  Memorex. 

Data  from  06C 

The  06C  subarray  is  again  connected  to  the  NORSAR  online  system. 
Unfortunately,  the  line  quality  is  not  good  enough  to  use  the  subarray 
as  a  part  of  the  online  processing.  Work  is  being  done  to  improve  the 
line  so  the  data  can  be  used  for  processing  purposes. 

Experience  with  bulletin  transfer 

The  new  bulletin  transfer  has  been  running  for  about  8  months,  and  our 
experience  is  very  good.  Even  with  transfer  3  times  a  day  our  work  is 
reduced  to  a  minimum  because  the  whole  process  is  automatic.  The 
transfer  in  weekend  and  holidays,  which  was  difficult  with  the  old 
system,  is  now  running  very  well. 

Future  plans  for  bulletin  transfer  are  to  use  ARPANET.  The  com¬ 
munication  between  the  USA  and  NORSAR  will  be  more  flexible.  Before 
ARPANET  is  used  for  bulletin  transfer,  some  problems  have  to  be 
cleared  up.  Problem  areas  are: 

-  regular  use  of  an  experimental  network  (the  SATNET  connection 
between  Europe  and  the  USA) 

availability  of  the  Norwegian  connection  to  SATNET. 

Today  we  have  an  ARPANET  connection  through  an  IBM/ PC.  The  IBM/ PC  is 
connected  to  Proteon's  PRONET  (10  Mbit  network)  which  again  is  inter¬ 
faced  to  SATNET. 


We  have  done  some  experimental  transfers  to  CSS.  A  new  ARPANET  inter¬ 
face  directly  connected  to  the  4341  is  preferred  for  a  regular 
transfer.  The  connection  will  eventually  be  made  through  a  Device 
Attachment  Control  Unit  (DACU)  from  IBM.  The  DACU  is  again  connected 
to  PRONET.  With  the  DACU  installed  NORSAR  can  again  be  a  full  host  on 
the  ARPANET. 


R.  Paulsen 


-  16  - 

V.  MAINTENANCE  ACTIVITIES 

V«1 _ Activities  in  the  Field  and  at  the  Maintenance  Center 

The  majority  of  the  visits  to  the  subarrays  in  the  period  have  been  in 
connection  with  corrective  work,  including  repair  of  cables.  Much  time 
has  been  spent  in  the  NORESS  array  covering  investigations  of  noise 
and  crosstalk  problems,  fitting  up  the  station,  clearing  up  and 
inspection  of  the  NORESS  area,  improving  the  local  oscillator  stabi¬ 
lity  and  reducing  noise  and  crosstalk. 

Otherwise,  repair  of  cards  and  other  electronic  equipment  has  been 
carried  out  at  the  NMC.  At  the  NDPC  different  SP-channel  parameters  in 
Che  NORSAR  array  have  been  evaluated,  after  data  collection  and 
analysis.  Remote  adjustments  of  the  mass  position  and  free  period  on 
LP-lnstruments  have  been  done. 

A  summary  of  the  activities  in  the  period  are  given  in  Table  V.l. 


Subarray/ 


Visit  to  the  subarray  in  connection  with 
errors  in  the  data  from  the  subarray 

Location  and  repair  of  SP02  cable 


Remote  Centering  Device  (RCD)  replaced  on 
EW  LP-seismometer.  Protection  card  SP04 
replaced 

Visit  to  the  subarray  in  connection  with 
deteriorated  data.  After  an  attempt  to  replace 
the  modem,  the  comm,  cable  was  found  to  be  the 
source  of  the  errors.  Repair  was  an  NTA 
responsibility 


3  Oct  84 


17,  24, 
25  Oct 

18  Oct 


18  Oct 


VVJSVJ 


Task 


Subarray/ 

area 


Date 


06C 

NORSAR 

array 

01A 

01B 

06C 

NORESS 

area 


NORSAR 

array 

01A 

NORESS 


NMC 

NORSAR 

array 


After  the  NORESS  activities  in  the  area,  15  Oct 
cables  had  to  be  spliced  in  the  period 

Under  remote  control  from  NDPC  17  mass  Oct 

position  and  2  free  period  adjustments  were 
done  on  totally  19  LP-seismometers  (of  which 
2  had  their  FP  adj.) 

SP  cable  (02)  repaired  after  breakage  6,7  Nov 

Located  damaged  SP  cable  (04).  Replaced  Remote 
Centering  Device  (RCD)  for  FP  and  Mass  Position 
(MP)  NS  seism. ,  R-A5  SP02,  05  adjusted 


Check-out  of  cables  and  electronic  equipment  12,  15 
in  connection  with  restart  of  subarray  Nov 

operation 

Final  inspection/clearing  of  the  area  after  16,  19-21 
the  field  and  construction  work.  Besides  23,  26-28 

engagement  in  connection  with  erecting  a  Nov 

fence  and  road  bar 

Monitoring  and  control  of  the  NORSAR  array  Nov 

(from  NDPC)  resulted  in  9  adjustments  of  MP 
and  2  FP  adjustments  on  totally  11 
LP-seismometers 


In  connection  with  private  trenching  work  in  5  Dec 

the  area  NMS  staff  had  to  point  out  NORSAR 

cables 

Quite  a  few  days  spent  on  different  tasks,  Dec 

such  as  fitting  up  the  central  station, 
clearing  the  station  area,  etc. 

Repaired  faulty  circuit  cards  on  various  Dec 

electronic  equipment 

Under  remote  control  from  NDPC  22  adjustments  Dec 
of  Mass  Position  (on  21  LP-seismometers),  and 
4  adjustments  of  Free  Period  (on  1  LP- 
seismometer)  were  carried  out 


I 


Task 


Subarray/ 

area 


Date 


02B 

Telemetry 

Battery  replacement  on  receiver  station 

2  Jan 

01A 

Replaced  RA-5  ampl.  on  SP01  and  00 

4  Jan 

01A 

Attenuated  channel  (04)  moved  from  00  to  the 

60  m  borehole 

OlA 

NMC  staff  to  the  subarray  in  connection  with 
line  problems 

8  Jan 

03C 

On  NS-seismometer  replaced  data  coil,  MP 
and  FP  motors.  Battery  charger  repaired,  dis¬ 
tilled  water  added  to  batteries.  Dry-out  of 
the  LPV  started  due  to  excessive  moisture 

17  Jan 

02B 

Telemetry 

Attempts  to  access  the  subarray  failed  due  to 
snow  problems 

29  Jan 

03C 

Inspection  of  the  LPV  still  indicated  moisture 

30  Jan 

06C 

New  NOKIA  modems  2400  baud  installed  in  the 

CTV  and  at  the  NDPC.  Also  PCM  line  reconfigured 
with  respect  to  equalizers.  Good  test  results 
modem-to-modem,  but  with  data  equipment 
(Modcomp/SLEM)  hooked  up,  error  at  an  average 
of  0.2  X  a  day  was  observed.  Investigation 
to  be  continued. 

30  Jan 

NORESS 

Due  to  noise  problems,  RC-f liters  were 
installed 

Jan 

NORSAR 

Array 

Remote  adjustment  (8)  of  Mass  Pos.  (on  6  LP 
seismometers).  In  addition  regular  checks  of 
Mass  Pos .  (MP)  and  Free  Period  (FP) 

Jan 

01B,  02B, 
02C,  03C 

Data  collection  and  analysis  in  connection 
with  discovering  possible  missing  numbers 
in  ADC  conversion 

Jan 

03C 

Dry-out  of  LPV  continued,  conclusion  after 
subarray  visit 

18  Feb 

NORESS 

NMC  staff  mainly  engaged  in  solving  noise  and 
cross-talk  problems 

Feb 

NORSAR 

array 

Remote  adjustment  (5)  of  Mass  Position  (5  LP 
seismometers  and  1  Free  Period  adjustment 

Feb 

-  19  - 


Subarray/ 

area 

Task 

Date 

01B,  04C 

Channel  parameters  as  RA-5  gain,  seismometer 
sensitivity  and  channel  resolution  calculated 
(by  means  of  the  CHANEV  SP  program) 

01A,  04C 

Data  collection  and  analysis  in  connection  with 
discovering  possible  missing  numbers  in  the  ADC 
conversion 

02B 

(Telemetry) 

NTA  Hamar  finished  the  work  on  the  telephone 
line  and  stations  back  in  operation 

12  Mar 

01A 

Subarray  visited  in  connection  with  SP(04) 
data  180°  out-of-phase 

19  Mar 

01A 

NMC  staff  to  the  subarray  to  point  out  an  SP 
cable  in  connection  with  trenching  work 

20  Mar 

NORESS 

Accuracy  of  the  Local  Oscillator  improved 
by  replacing  the  LO  crystal 

March 

NORESS 

Station  D4  visited  and  Fiberoptic  Transmitter 
replaced 

March 

NORSAR 

array 

Remote  adjustment  (10)  of  Mass  Pos.  (on  8  LP 
seismometers)  and  Free  Period  (3)  on  same 
number  of  LP  seismometers. 

March 

Table  V.l 


Activities  in  the  field  and  at  the  NORSAR  Maintenance 
Center,  October  1984  -  March  1985. 


V.2 _ Improvements  and  Modifications 

Attenuated  channel  (04)  moved  from  borehole  00  to  the  60  m  hole. 

V.3 _ Array  Status 

As  of  31  March  1985  the  following  channels  deviated  from  tolerances 

01B  05 
04C  03 

06C  (not  in  operation  due  to  comm,  system) 

01A  01  8  Hz  filter 

02 

04  Attenuated  30  dB 


0. A.  Hansen 
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VII. 
VII. 1 


SUMMARY  OF  TECHNICAL  REPORTS/ PAPERS  PREPARED 

Real  time  event  detection  using  the  small-aperture  NORESS 


Array 

Since  January  1985  data  from  the  new  small-aperture  array  NORESS  In 
Norway  have  been  processed  In  real  time  at  the  NORSAR  data  center  at 
Kjeller.  The  data  used  In  the  detection  processing  comprise  25  SPZ 
channels,  deployed  over  an  area  3  km  in  aperture  and  sampled  at  a 
40  Hz  rate.  The  detection  algorithm  has  been  described  by  Mykkeltveit 
and  Bungum  (1984),  and  briefly  consists  of 

-  Digital  narrow-band  filtering  (six  filters) 

Beamforming  (conventional  and  incoherent) 

STA/LTA  detector  applied  to  each  beam 

Frequency-wavenumber  analysis  of  detected  signals 

Association  of  regional  phases  to  aid  in  locating  events. 

In  the  following,  initial  results  from  this  processing  are  described, 
together  with  recommendation  for  future  research. 

Teleseismic  events 

The  P-wave  signal-to-noise  ratio  at  NORESS  typically  peaks  in  the 
2-4  Hz  band  for  teleseismic  events  in  Eurasia,  whereas  lower  dominant 
frequency  is  often  observed  for  signals  from  the  western  hemisphere 
and  for  epicentral  distances  exceeding  70  degrees.  Signal  coherency 
is  naturally  very  good  for  teleseismic  P  across  NORESS,  even  at  very 
high  frequencies.  Depending  on  signal  frequency,  the  best  beam  SNR  is 
obtained  using  various  subsets  of  the  array;  e.g.,  at  2  Hz  the  best 
subset  consists  of  the  center  instrument  AO  together  with  the  C  and  D 
rings  (see  Fig.  VII. 1.1).  Using  this  subset,  the  beamforming  gain  meets 


or  exceeds  /N,  due  Co  Che  noise  suppression  characCerisCics  described 
in  SecCion  VII. 2. 

An  example  of  Celeseismic  processing  results  from  NORESS  is  given  in 
Fig.  VII. 1.2.  The  improvemenC  in  SNR  on  Che  beam  relaCive  Co  Che 
single  sensors  is  quiCe  remarkable. 

The  on-line  procedure  uses  sCeered  Celeseismic  beams  raCher  chan  a 
single  inf iniCe-velociCy  beam  for  each  filcer.  The  resulClng  SNR  gain 
is,  as  an  example,  4-5  dB  for  an  apparenC  velociCy  of  16  km/s,  aC  3  Hz 
frequency.  For  lower  frequencies  or  higher  velociCles,  Che  gain  is 
less.  Thus,  applying  only  inf inCe-velociCy  beamforming  for  Celeseismic 
dececClon  is  a  viable  alCernadve,  in  view  of  Che  possibiliCy  of 
operaCing  at  a  lower  detection  threshold  with  Che  same  false  alarm 
race. 

The  NORESS  array  detects  many  Celeseismic  signals  not  observed  at  Che 
large  aperture  NORSAR  array,  which,  like  NORESS,  is  located  in 
southeastern  Norway;  especially  from  selected  regions  in  Eurasia.  On 
Che  other  hand,  Che  signal  focusing  effects  underneath  NORSAR  cause 
some  instruments  to  have  up  to  an  order  of  magnitude  stronger  signals 
Chan  NORESS  sensors,  for  regions  such  as  Hindu  Kush  and  Che  Kuriles. 
For  these  regions,  NORESS  does  not  match  the  NORSAR  detection 
capability. 

Because  of  the  small  aperture  of  NORESS,  only  a  very  coarse  automatic 
estimate  of  the  location  of  teleselsmic  events  is  currently  being 
made.  Azimuth  errors  depend  on  phase  velocity  and  frequency,  and  are 
often  around  5-10  degrees  for  small  events.  Slowness  estimates  from 
the  automatic  process  are  typically  about  1  sec/deg  different  from 
those  of  NORSAR.  Nevertheless,  it  is  clear  that  location  estimates 
from  a  small  array  like  NORESS  will  be  very  valuable  to  provide  a 
starting  point  for  association  procedures  using  a  seismic  network. 


Regional  correcCions  and  more  detailed  off-line  analysis  are  required 
to  assess  the  eventual  capabilities  of  the  array  in  this  regard. 

Local  and  regional  events 

At  local  and  regional  distances,  the  best  SNR  for  the  P  wave  varies 
from  3-5  Hz  (at  around  1000  km)  to  more  than  8  Hz  (local  distances). 
Consequently,  either  steered  beams  or  incoherent  beamforming  is 
necessary  to  exploit  the  array  capability.  P-signal  coherency  at 
NORESS  is  usually  good  enough  to  utilize  the  full  array  for  F-K  pro¬ 
cessing  of  local  and  regional  P-phases,  at  least  up  to  about  6-8  Hz. 

The  Lg  phase  is  usually  of  slightly  lower  frequency  than  P.  Conven¬ 
tional  beamforming  is  not  very  efficient  for  Lg,  since  the  preceding 
coda  (from  Sn)  comes  in  with  about  the  same  phase  velocity  and 
azimuth.  Consequently,  little  "noise"  suppression  takes  place.  A  pro¬ 
mising  approach  is  that  of  performing  narrow-band  filtering  at  several 
frequencies  for  the  purpose  of  Lg  detection.  It  turns  out  that  the  SNR 
is  greatly  Improved  in  those  frequency  bands  where  P  and  Sn  coda 
energy  are  low  relative  to  the  Lg  energy.  Combining  narrow-band 
filtering  with  incoherent  beamforming  has  been  found  particularly 
effective. 

An  example  of  a  complete  record  from  a  regional  event  processed  at 
NORESS  is  given  in  Fig.  VII. 1.3.  The  detection  times  for  P  and  Lg  are 
marked  on  the  panel  of  Fig.  VII. 1.3a,  whereas  Fig.  VII. 1.3b  shows  F-K 
solutions  for  P  and  Lg  together  with  short  plots  of  each  phase  in  an 
expanded  time  scale.  The  estimated  azimuths  for  P  and  Lg  differ  by  2 
degrees  in  this  case;  and  a  deviation  of  0-5  degrees  is  common. 
However,  a  difference  of  10  degrees  and  more  is  also  fairly  often 
observed  in  the  automatic  solution. 
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The  location  accuracy  of  NORESS  for  regional  events  is  currently  being 
studied.  No  statistically  reliable  results  are  available  so  far  due  to 
the  limited  data  base  of  known  locations.  The  procedure  of  locating 
events  by  associating  P  and  Lg  is  applicable  only  up  to  about  1200  km 
distance.  At  greater  distances,  the  Lg  usually  is  too  weak  to  be 
detected  automatically,  but  can  sometimes  be  identified  by  visually 
inspecting  the  waveform  plots. 

In  conclusion,  the  initial  results  obtained  from  NORESS  are  very  en¬ 
couraging,  and  have  met  and  in  some  cases  exceeded  the  expectations. 
Data  quality  has  been  excellent,  and  the  operational  stability  has 
been  very  satisfactory,  taken  into  account  the  complexity  of  the 
system.  Particularly  noteworthy  is  the  excellent  P  detection  in  the 
2-4  Hz  band,  which  is  due  to  greater  than  /N  noise  suppression  com¬ 
bined  with  strong  P-wave  energy.  The  automatic  detection  of  secondary 
phases  needs  refinement,  and  in  particular  the  narrow-band  filter  bank 
processing  should  be  further  investigated. 


F.  Ringdal 
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Fig.  VII. 1.2  Example  of  automatic  processing  of  a  teleseismic  event 
using  NORESS.  The  on-line  F-K  solution  is  shown  at  the 
top,  together  with  detection  parameters.  Ten  individual 
sensor  traces  (filtered  2. 5-4. 5  Hz)  are  plotted, 
together  with  the  array  beam  (bottom). 
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Fig.  VII. 1.3b  NORESS  automatic  analysis  results  for  the  event  of  Fig. 

VII. 1.3a.  The  P-phase  (left)  and  Lg  phase  (right)  are 
shown  for  all  SPZ  instruments.  F-K.  solutions  are  shown 
for  P  (top)  and  Lg  (bottom). 
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VII. 2 _ Preliminary  evaluation  of  the  event  detection  and  location 

capability  of  the  small-aperture  NORESS  array 
The  new  NORESS  array  In  Norway  has  been  operation  since  October  1984. 
The  geometry  of  the  array  Is  shown  In  Fig.  VII.2.1.  The  array  has  25 
short  period  vertical  seismometers,  arranged  on  concentric  rings 
around  the  central  site.  Four  of  the  25  sites  are  occupied  by  three- 
component  seismometers. 

A  data  base  of  local  and  regional  events  has  been  collected,  with  the 
purpose  of  evaluating  the  event  detection  and  location  capability  of 
the  new  array.  As  a  first  approach,  we  have  concentrated  our  efforts 
on  the  study  of  a  data  set  of  18  seismic  events  that  occurred  in  the 
Baltic  and  Leningrad  regions  of  the  U.S.S.R.  during  the  time  period 
January-March  1985.  The  events  are  associated  with  mining  activity  and 
are  all  reported  and  located  by  the  Finnish  network  of  seismic  sta¬ 
tions.  An  example  of  one  such  event  is  shown  in  Fig.  VII. 2. 2.  Three 
distinct  phases  P,  Sn  and  Lg  can  be  seen  for  all  18  events.  Spectra 
for  noise  preceding  the  P  onset,  the  P  phase,  the  Sn  phase  and  the  Lg 
phase  are  plotted  together  in  Fig.  VII. 2. 3  for  one  of  the  events.  As 
can  be  seen,  the  P  phase  is  very  rich  in  high  frequencies,  while  Sn 
and  Lg  approach  the  noise  level  from  about  8  Hz.  Note  also  the  noise 
peak  at  around  6  Hz,  which  is  typical  of  daytime  noise  samples. 

Now,  to  assess  the  detection  capability  of  the  array,  we  have  esti¬ 
mated  the  SNR  gain  by  beamforraing  for  the  three  phases  P,  Sn  and  Lg. 
This  has  been  done  by  computing  the  signal  loss  and  noise  suppression 
and  taking  the  ratio.  The  signal  loss  is  the  spectrum  of  the  beam 
(with  appropriate  bearing  and  velocity  depending  on  the  phase  type) 
divided  by  the  average  spectrum  of  all  individual  channels.  Similarly, 
the  noise  suppression  is  the  ratio  of  the  spectrum  of  the  noise  beam 
(noise  preceding  the  signal,  shifted  in  the  same  way  as  the  signal 
itself  for  the  beamforming)  to  the  average  spectrum  of  all  individual 
channels. 
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The  noise  spectra  have  been  estimated  using  the  indirect  covariance 
method.  We  first  estimate  the  correlation  function  by  splitting  a  long 
data  record  into  many  windows,  calculating  a  sample  correlation  func¬ 
tion  for  each  window,  then  averaging  the  sample  correlation  functions. 
Typically,  we  use  20  windows,  each  of  which  is  5  seconds  long.  Because 
the  earth  noise  has  such  a  large  dynamic  range,  we  prewhiten  it  prior 
to  estimating  the  correlation  function  with  a  low-order  prediction- 
error  filter.  The  spectrum  is  then  estimated  by  windowing  the  correla¬ 
tion  function  with  a  3  second  Hamming  window,  then  computing  the 
Fourier  Transform.  The  spectral  estimate  obtained  this  way  is  compen¬ 
sated  then  for  the  effects  of  prewhitening. 

The  signal  spectrum  is  calculated  by  windowing  the  signal  with  a  10Z 
cosine  taper  window,  5  seconds  long,  then  calculating  the  Fourier 
Transform  and  finally  squaring  the  transform. 

When  we  do  direct  comparison  of  the  noise  spectrum,  which  is  a  power 
density  spectrum,  with  the  signal,  which  is  an  energy  density 
spectrum,  we  divide  the  signal  spectrum  by  the  length  of  the  signal 
analysis  window  to  convert  energy  density  to  power  density.  The  two 
spectral  quantities  are  then  directly  comparable.  This  has  been  done 
in  Fig.  VII. 2. 3. 

The  signal  loss,  noise  suppression  and  SNR  gain  by  beamforming  have 
been  computed  for  four  different  array  configurations: 

All  elements  Included  ("ALLV"  in  the  figures) 

All  elements  minus  A-ring  ("INTERMEDIATE") 

Center  instrument,  C-ring  and  D-ring  ("TELEV") 

-  All  elements  minus  D-ring  ("HIFREQ"). 

In  Fig.  VII. 2. 4,  noise  suppression  curves,  representing  averaging  over 
the  18  events,  are  shown  for  all  four  array  geometries.  The  horizontal 
lines  represent  /N  suppression  of  noise,  where  N  is  the  number  of  sen¬ 
sors  in  the  geometry  under  consideration.  This  is  the  expected  level 


for  noise  that  is  uncorrelated  over  the  entire  array.  As  can  be  seen, 
the  "TELEV" -geometry  is  remarkably  good  in  suppressing  the  noise  in 
the  1.5-3. 5  Hz  band,  while  the  other  geometries  approach  the  /N  level 
for  various  frequencies,  but  do  not  show  better  noise  suppression  than 
/N.  These  results  partly  reflect  on  the  SNR  gains  shown  in  Fig. 

VII. 2. 5.  the  signal  losses  (not  shown)  tend  to  be  larger  for  the 
larger  geometries  (including  outer  rings)  because  of  the  decaying 
signal  correl  tion  with  increasing  station  separation.  The  overall 
effect  is  to  level  the  gain  curves,  and  we  can  see  from  Fig.  VII. 2. 5 
that  only  at  fairly  low  and  fairly  high  frequencies  do  the  gain  curves 
for  the  different  geometries  deviate  substantially. 

For  the  Sn  and  Lg  phases,  the  "noise  suppression"  estimates  reflect 
suppression  by  beamforming  of  the  P  and  Sn  coda,  respectively,  i.e., 
we  address  the  problem  of  detecting  a  secondary  phase  in  the  coda  of  a 
preceding  phase.  Results  shown  in  Fig.  VII. 2. 6  are  for  the  Lg  phase  and 
correspond  to  the  best  SNR  result  obtained  among  the  four  geometries. 
SNR-gain  is  now  well  below  /N,  essentially  because  of  failure  to 
suppress  the  Sn  coda,  which  propagates  with  similar  phase  velocity  and 
azimuth  as  the  Lg  onset.  Still,  there  is  a  6  dB  gain  from  beamforming 
at  1.2  Hz. 

Fig.  VII. 2. 7  shows  the  P  onset  for  a  Novaya  Zemlja  explosion  (distance 
approx.  20°)  recorded  on  NORESS,  and  the  SNR  gain  by  beamforming  com¬ 
puted  as  before  for  the  "TELEV"  configuration.  In  the  frequency  range 
1.0-2. 7  Hz  the  gain  is  better  than  /N,  because  of  a  very  modest  signal 
loss  in  combination  with  optimum  noise  rejection,  as  for  the 
"TELEV"-conf iguration  in  Fig.  VII. 2. 3.  The  same  result  is  visualized 
in  the  time  domain  in  Fig.  VII. 2. 8. 

Our  approach  to  the  location  performance  evaluation  has  been  to  com¬ 
pute  narrow-band  FK-spectra  for  the  P,  Sn  and  Lg  signals  for  the 
events  in  the  data  base.  FK-spectra  have  been  computed  for  the  array 
subsets  given  above  and  for  frequencies  corresponding  to  distinct 
spectral  maxima  for  each  phase.  Typical  results  are  given  in 
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Table  VII. 2.1  for  one  event.  As  can  be  seen,  the  analysis  frequency  is 
the  critical  parameter,  and  more  so  than  the  array  subgeometry  used  in 
the  evaluation  of  the  FK-spectra.  According  to  these  results,  then, 
the  arrival  azimuth  varies  fairly  strongly  with  frequency  and  due  care 
must  be  taken  in  selecting  the  analysis  frequency.  To  gain  more 
insight  into  these  frequency-dependent  lateral  refraction  effects,  we 
have  computed  the  broad-band  maximum  likelihood  FK-spectra  for  a 
number  of  phases.  One  example  is  shown  in  Fig.  VII. 2. 9,  where  the  azi¬ 
muth's  dependency  on  frequency  is  given  implicitly  in  the  contour 
plot.  The  azimuth  value  of  89.5  serves  to  illustrate  that  the  broad¬ 
band  FK  tends  to  give  more  stable  estimates  than  those  derived  from 
narrow-band  FK-spectra. 

The  research  effort  continues  with  analysis  of  events  from  other 
regions,  in  order  to  obtain  the  capabilities  of  the  new  NORESS  array 
for  local  and  regional  events  from  all  source  regions  of  Interest.  As 
experience  is  gained,  the  results  from  this  study  will  be  utilized 
directly  in  the  online  processing  of  NORESS  data. 


S.  Mykkeltveit 

D.B.  Harris,  Lawrence 

Livermore  National  Laboratory 

T.  Kvaerna 
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VII. 2.5  SNR  gain  by  beamforming  of  regional  F  phase  for  the  four 

array  subgeometries.  The  horizontal  line  represents  VN  gain, 
where  N  Is  the  number  of  sensors  in  each  subgeometry.  The 
vertical  scale  is  logarithmic. 
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Fig.  VII. 2. 6  Noise  suppression  and  SNR  gain  for  the  Lg  phase  for  sub¬ 
geometry  "TELEV".  Horizontal  lines  represent  /N-noise 
rejection  (upper  plot)  and  /N-beamforming  gain  (lower  plot), 
with  N  ■  17  for  this  subgeometry.  Vertical  scales  are 
logarithmic. 
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Fig.  VII. 2. 8  Time  domain  illustration  of  the  SNR-gain  results  in  Fig. 

VII. 2. 7.  Top:  Normalized  plot  of  three  single  channels 
(out  of  the  17  participating  in  the  beam)  and  the  beam, 
showing  a  very  modest  signal  loss.  Bottom:  Noise  rejec¬ 
tion  obtained  on  preceding  noise  (true  amplitude  plot). 
The  filter  bandpass  is  1.3-2. 5  Hz,  corresponding  to  the 
peak  in  the  gain  curve  in  Fig.  VI I. 2. 7. 
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Fig.  VII. 2. 9  Broad-band  FK-analysis  for  the  Lg  phase  of  the  event  used 
in  Table  VII. 2. I.  In  order  to  obtain  the  velocity  corre¬ 
sponding  to  the  peak,  the  value  of  3.15  km/s  must  be 
multiplied  by  the  dominant  frequency  of  the  Lg  signal 
(1.33  Hz)  to  yield  a  phase  velocity  of  4.19  km/s. 
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VII. 3  Seismic  modelling  of  an  anisotropic  boundary  layer 
Various  parts  of  the  upper  mantle  have  been  proposed  to  be 
anisotropic.  More  recently  it  has  been  suggested  that  anisotropy  may 
also  be  present  at  other  depths,  especially  in  boundary  layers 
experiencing  differential  flow.  D"  is  such  a  layer.  More  specifically, 
we  may  need  to  consider  a  model  of  D"  with  two  distinct  structures, 
resembling  the  continent-ocean  division  of  the  lithosphere.  The  possi¬ 
bility  of  observing  the  appearance  of  a  layered  structure  would  be 
restricted  to  the  permanent  (continental-type)  zones,  whereas  intrin¬ 
sic  anisotropy  predicted  by  boundary  layer  theory  would  prevail  in  the 
transient  (ocean-floor  type)  zones.  The  possibility  of  lateral  hetero¬ 
geneity  has  an  added  significance  in  this  context  because  it  has  been 
suggested  that,  when  interpreted  in  terms  of  one-dimensional  models, 
such  structures  may  appear  anisotropic  (Crampin  et  al,  1984).  For  the 
purpose  of  seismic  modelling,  we  assume  uniaxial  anisotropy  (with  sym¬ 
metry  axis  vertical),  often  termed  transverse  isotropy  in  the  seismo- 
logical  literature,  not  that  we  believe  this  to  apply  locally,  but 
because  we  assume  that  azimuthal  anisotropy  is  averaged  out  by  the 
several  wave  paths  used  in  observational  work.  The  required  modifica¬ 
tions  to  reflectivity  type  of  methods  have  been  discussed,  a.o.  in  a 
recent  contribution  (Doornbos  et  al,  1985).  Piecewise  smooth  models 
were  considered  here.  The  formal  expressions  for  (generalized)  reflec¬ 
tion  coefficients  are  unchanged,  but  the  required  modifications  con¬ 
cern  the  fundamental  matrices,  and  the  vertical  slownesses  and  their 
integrals  over  radius,  the  so-called  "Tau  functions”. 

One  Inference  is  that  for  radius  r  not  far  from  the  turning  point  rQ, 
x  is  determined  mainly  by  the  profile  of  horizontal  wave  velocity. 

This  is  in  accordance  with  the  results  of  a  numerical  experiment 
giving  the  effect  of  anisotropic  D”  on  the  diffracted  wave  fields  of  P 
and  SH.  In  Fig.  VII. 3.1  are  shown  the  changes  (relative  to  the  isotro¬ 
pic  approximation)  in  logarithmic  attenuation  with  distance  and  dT/dA, 
as  functions  of  frequency,  which  is  the  form  in  which  most  obser¬ 
vational  data  have  been  presented.  From  these  results  we  conclude  that 
diffracted  P  and  SH  are  mainly  controlled  by  the  horizontal  velocities 


ol  and  pL,  respectively,  and  these  are  related  to  the  elastic 
constants  A  and  N  of  a  uniaxial  structure.  On  the  other  hand,  five 
elastic  constants  (A,C,N,L  and  F)  are  needed  to  fully  characterize 
this  structure,  and  four  constants  (A,C,N  and  F)  are  needed  to  obtain 
the  bulk  modulus  K.  What  is  usually  inferred  for  the  purpose  of  tem¬ 
perature  calculations  is  an  apparent  bulk  modulus  K'  based  on  the 
assumption  of  an  isotropic  structure:  X.+2p  «*  A,  ^  ■  N.  We  find 

3  1-N/A-(F/A)2  1-C/A 

K/K.  „  -  .  - - (1) 

(3-4N/A)  2-N/A-2F/A  3(3-4N/A) 


where 

F/A 


N 

nd-2  - 
A 


L 

~) 

N 


Results  based  on  eq.  (1)  are  given  in  Table  V11.3.I.  These  results  can 
be  used  to  conclude  that  a  small  amount  of  anisotropy  requires  a 
correction  to  the  selsmlcally  inferred  temperature  increment  based  on 
K' .  For  D"  it  is  estimated  that  the  correction  can  be  up  to  about 


400°K. 


D.J.  Doornbos 
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Table  VII. 3.1  Ratio  of  uniaxial  and  isotropic  bulk  modulus  K/K' 

Assumed  relation  between  Isotropic  and  uniaxial  elastic 
constants: 

A  +  2)i  -  A,  p  ■  N 

Results  based  on  eq.  (1)  with  N/A  ■  0.28  and  C/A  *  1. 

If  C/A  *  1,  subtract  (l-C/A)/5.64  from  all  results. 

The  broken  line  encloses  the  ranges  of  elastic 
constants  that  appear  to  us  to  be  the  most  plausible. 
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Fig.  VII. 3.1 


1/64  1/16  1/8 
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Effect  of  a  transversely  isotropic  layer  on  logarithmic 
attenuation  (y)  and  dT/dA  of  diffracted  P  and  SH  around 
the  core.  The  distance  range  considered  is  100-135°. 

The  velocities  a  and  p  in  the  isotropic  part  of  the 
model,  and  the  horizontal  velocities  for  P  and  p^ 
for  SH  in  the  layer,  correspond  to  model  A  of  Doornbos 
et  al  (1985).  The  anisotropy  is  zero  at  the  top  of  the 
layer,  and  the  maximum  at  the  bottom.  The  maximum  ani¬ 
sotropy  for  SH  is  Ph/Pv  "  1.02.  The  three  P-wave  ani¬ 

sotropies  considered  are  (1)  a^/a^  -  1.01,  ti  *  1, 

Ph/Pv  -  1;  (2)  och/av  -  1,  r,  -  0.95,  ph/pv  -  1;  (3) 

ach/av  ■  1,  tj  ■  1,  Ph/Pv  ■  1.02.  The  values  plotted  are 

differences  from  results  with  isotropic  model  A. 


VII. 4  An  application  of  long-period  waveform  inversion  to  shallow 

sources  under  the  Tibetan  Plateau 

In  previous  semiannual  summaries  we  have  reported  on  the  development 

of  long-period  waveform  inversion  methods.  We  have  now  applied  these 

methods  to  analyze  shallow  focus  sources  under  the  Tibetan  Plateau. 

The  purpose  of  this  work  was  threefold: 

1)  To  compare  the  performance  of  different  methods  for  deter¬ 
mining  source  parameters  and  location,  and  for  determining 
travel  time  and  amplitude  residuals 

2)  To  test  current  hypotheses  concerning  recent  tectonics  of 
Tibet,  as  implied  by  the  principal  stress  pattern  and  the 
source  depths 

3)  To  determine  station  residuals  and  compare  these  residuals 
to  those  for  deep  events. 


The  methods  we  have  applied  include  a  centroid-moment  tensor  method 
comparable  to  the  one  described  by  Dzlewonski  and  Woodhouse  (1983), 
and  a  correlation  method  as  described  by  Doornbos  (1983).  Both  methods 
require  the  calculation  of  Green's  functions  in  a  realistic  earth 
model;  we  have  used  a  WKBJ  method  for  mixed  paths,  and  we  have 
modified  PREM  by  Incorporating  in  the  source  region  a  crustal  model 
for  Tibet,  and  under  the  receivers  an  appropriate  continental  model. 
The  principal  differences  between  the  methods  concern  the  nature  of 
the  inversion  procedure  (linearized  versus  non-linear  inversion)  and 
the  optimization  criterion  (least-squares  error  versus  maximum 
correlation).  On  the  basis  of  analysis  of  GDSN  data  with  P  and  SH 
waveforms  from  seven  events  in  Tibet,  we  conclude  that  the  non-linear 
correlation  method  leads  to  more  stable  results.  There  are  several 
reasons  for  this:  (1)  the  correlation  method  is  less  sensitive  to  the 
Initial  value  for  source  location.  This  is  an  important  point  since 
the  standard  focal  depth  value  (33  km)  used  by  NE1S  for  most  of  these 
events  is  significantly  more  than  the  values  obtained  in  this  work 


(<10  km).  (2)  The  correlation  method  explicitly  accounts  for  travel 
time  anomalies  and  is  robust  with  respect  to  amplitude  anomalies  at 
the  different  stations.  An  example  of  a  result  by  the  correlation 
method  is  given  in  Fig.  VII.4.1.  The  fault  plane  solution  in  this 
figure  illustrates  that  faulting  is  consistent  with  east-west  exten- 
sional  tectonics.  This  was  in  fact  found  for  all  events  as  summarized 
in  Fig.  VII. 4. 2,  and  it  supports  one  of  the  competing  theories  which 
have  been  proposed  for  this  region.  The  fact  that  all  focal  depth 
values  obtained  by  the  correlation  method  are  less  than  10  km  is 
moreover  consistent  with  a  relatively  hot  crust  (c.f.  Chen  and  Molnar, 
1983).  In  the  previously  reported  analysis  of  deep  events  it  was  found 
that  the  travel  time  residuals  of  S  are  surprisingly  large;  there  was 
both  a  base  line  effect  and  a  trend  with  epicentral  distance.  It  is  of 
interest  to  compare  this  result  to  the  travel  time  residuals  for  the 
shallow  focus  events  analyzed  here.  These  residuals  are  plotted  in 
Fig.  VII. 4. 3;  a  correction  for  finite  source  duration  has  been  applied 
here.  To  interpret  the  P  residuals,  note  that  the  origin  times  given 
by  NEIS  involve  the  J-B  velocity  model  and  focal  depths  which  are 
usually  taken  to  be  33  km.  Travel  time  effects  due  to  systematic  error 
in  focal  depth  and  due  to  the  difference  between  the  J-B  model  and  the 
modified  PREM  model  used  here  amount  to  about  -3.5  s.  This  explains 
the  P  residuals  in  Fig.  VII. 4. 3,  as  expected. 

The  S  residuals  would  have  to  be  correspondingly  corrected  by  adding  + 
3.5  8.  The  corrected  S  residuals  at  the  larger  epicentral  distances 
are  then  slightly  negative,  but  on  average  the  residuals  are  small. 
This  contrasts  with  the  large  positive  S  residuals  for  deep  events.  It 
suggests  the  need  for  modifications  of  the  PREM  model  involving  both 
the  upper  and  lower  mantle. 

0.  Pettersen 

D.J.  Doornbos 
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Fig.  VII. 4.1  Source  analysis  for  one  of  the  events  in  this  study. 

Observed  and  synthetic  waveforms  at  GDSN  stations, 
together  with  the  double-couple  approximation  to  the 
moment  tensor  solution,  in  equal  area  projection  of  the 
lower  focal  hemisphere.  The  centroid  depth  is  5  km.  The 
time  axis  is  in  10  s  intervals. 


A  =T  D=P  0=B 


Fig.  VII. 4. 2  Principal  axis  of  7  shallow  sources  under  the  Tibetan 
Plateau,  in  equal  area  projection  of  the  lower  focal 
hemisphere. 
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Fig.  VII. 4. 3  Travel  time  residuals  versus  distance  for  P  (a)  and  SH 
(b).  The  residuals  are  defined  as  observed  minus  pre¬ 
dicted.  Each  station  is  positioned  at  the  average 
distance  from  the  events  which  were  recorded. 


VII. 5 _ Interactive  analysis  program  for  use  with  NORESS  data 

We  have  completed  a  flexible  analysis  program  for  processing  the 
NORESS  data  in  an  Interactive  manner.  This  code  will  be  used  in  coming 
months  to  assess  the  performance  of  the  new  regional  array  and  to  exa¬ 
mine  a  variety  of  aspects  of  wave  propagation  in  the  region. 

The  program  provides  basic  file  manipulation  and  plotting  functions. 
Data  can  be  read  into  the  program  from  standard  NORESS  tapes  or  from 
disk  file  images  of  the  tapes.  The  program  also  has  the  capability  of 
writing  and  reading  data  to  disk  files  in  its  own  format  for  temporary 
storage.  The  data  storage  part  of  the  program  is  organized  as  a  signal 
stack,  that  is,  it  is  a  first-in  first-out  push  down  stack,  where 
entries  on  the  stack  are  signals.  Signals  may  be  popped  off  of  the  top 
of  the  stack  into  the  local-format  files,  or  pushed  onto  the  stack 
from  those  files. 

The  top  of  the  stack  is  the  reference  point  for  all  commands.  The 
default  operand  for  all  unary  operations  is  the  top  signal  on  the 
stack,  and  for  binary  operations  the  default  operands  are  the  top  two 
signals.  Usually,  a  command  can  be  extended  to  apply  to  any  number  of 
operations,  such  as  scaling,  absolute  value,  square  root,  are  provided 
as  well  as  binary  operations  like  addition,  subtraction,  multiplica¬ 
tion,  division  of  signals.  These  are  polntwise  operations.  Basic 
filtering,  spectral  analysis,  and  array  processing  functions,  such  as 
stacking  and  FK  calculations,  are  also  incorporated. 

Several  canned  plotting  routines  are  available  to  plot  trace  data, 
array  response  patterns  and  FK  spectra  as  contour  maps  or  perspective 
drawings,  power  spectra,  and  maps  of  array  geometry  (Fig.  VII. 5.1  and 
VII. 5. 2).  To  a  limited  extent,  the  plots  can  be  titled  and  labeled. 

The  viewing  window  for  the  trace  data  plotting  command  is  under  cursor 
as  well  as  keyboard  control.  This  provides  for  convenient  zooming  on  a 


particular  window  and  also  for  cursor  selection  of  the  analysis  window 
for  analysis  functions  that  operate  on  some  restricted  domain  of 
interest. 

The  code  is  command  driven,  rather  than  menu  driven,  to  enhance  the 
speed  of  interaction.  A  help  command  has  been  provided  to  list  the 
available  command  options  and  syntax.  A  primitive  macro  facility 
exists  as  well,  permitting  groups  of  commands  to  be  stored  in  external 
files  and  accessed  with  a  single  macro  call  (Fig.  VII. 5. 3).  Macro 
calls  can  be  nested  to  a  maximum  depth  of  four,  but  may  not  be 
recursive. 

A  facility  for  storing  temporary  numerical  and  character  values  has 
been  Implemented.  We  have  called  this  facility  the  "blackboard", 
although  it  is  far  from  a  blackboard  in  the  expert  system  sense  of  the 
term.  Blackboard  entries  are  named,  and  they  may  be  referred  to  by 
name  as  arguments  to  all  of  the  commands  in  the  program.  Simple  arith¬ 
metic  and  string  manipulations  can  be  performed  on  blackboard  entries 
as  well. 

The  purpose  of  the  blackboard  is  to  facilitate  automatic  processing  of 
large  amounts  of  data.  Blackboard  entries  are  useful  for  passing  argu¬ 
ments  to  macros,  effectively  turning  the  macros  into  subroutines.  This 
approach  makes  it  possible  to  build  up  extensive  and  fairly  sophisti¬ 
cated  signal  processing  and  measurement  capabilities  from  the  primi¬ 
tive  operations  supplied  in  the  program.  Our  Intent  was  to  provide  a 
flexible  "breadboard"  for  trying  new  ideas  that  would  also  have  the 
capability  for  production  work  with  a  new  method,  once  that  method  had 
proved  useful. 

D.B.  Harris,  Lawrence 
Livermore  Ntl.  Laboratory 

T.  Kvmrna 
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MAP  OF  ARRAY  ELEMENT  LOCATIONS 


1  KILOMETER 


Fig.  VII. 5.1  a)  Locations  of  the  vertical  instruments  of  the  NORESS 


EORETICAL  BEAM  PATTERN 

YEN  IN  DECIBEL  BELOW  MAX.  PEAK 


beam  pattern  of  the  NORESS  array. 

EORETICAL  BEAM  PATTERN 

YEN  IN  DECIBEL  BELOW  MAX.  PEAK 


Lew  of  the  theoretical  beam  pattern 
array. 
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Fig.  VII. S. 2  b)  F-k  solution  of  the  Lg  phase  of  the  event  in  Fig 
VII. 5. 2  a). 
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"  MACRO  TO  PREWHITEN  TOP  STACK  SIGNAL  AND  CALCULATE  ITS  SPECTRUM  WITH 
M  THE  EFFECTS  OF  PREWHITENING  COMPENSATED  FOR.  EFFECTIVELY 

"  FLUSHES  THE  TOP  STACK  SIGNAL.  SPECTRUM  ON  TOP  OF  STACK. 

1 

2  PREWHITEN  . 

3  SWAP  . 

4  SETW  ON  : SIGNAL-START  : SIGNAL-STOP  . 

5  CUT  . 

6  SETW  OFF  . 

7  WINDOW  COSINE  . 

8  PAD  TO  1024  . 

9  FFT  . 

10  SQR  2  . 

11  ADD  . 

12  SWAP  . 

13  FLUSH  . 

14  SWAP  . 

15  FFT  . 

16  SQR  2  . 

17  ADD  . 

18  SWAP  . 

19  FLUSH  . 

20  SWAP  . 

21  DIVIDE  . 

22  TERMINAL  . 


LOGIC  OF  THE  PW-SPECT  MACRO 


1  SIGNAL  ON  TOP  OF  THE  STACK 
STACK:  SIGNAL 


Fig.  VII. 5. 3  PW  spectra  (page  1  of  3) 
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2  PREWHITEN  . 

PREWHITEN  THE  THE  TOP  SIGNAL  AND  PUT  THE  PREWHITENING  FILTER  ON 
TOP  OF  THE  STACK  (1024  POINTS). 

STACK:  PREWHITENING  FILTER,  SIGNAL 

3  SWAP  . 

SWAP  THE  TWO  FIRST  SIGNALS 
STACK:  SIGNAL,  PREWHITENING  FILTER 

4  SETW  ON  : SIGNAL-START  : SIGNAL-STOP  . 

SET  ANALYSIS  WINDOW  ON  THE  SIGNAL  TRACE  .  THE  VALUES  OF  SIGNAL-START 
AND  SIGNAL-STOP  ARE  CALLED  FROM  THE  BLACKBOARD  . 

STACK:  SIGNAL,  PREWHITENING  FILTER 

5  CUT  . 

STORE  THE  PART  OF  THE  SIGNAL  WITHIN  THE  ANALYSIS  WINDOW,  ERASE  THE 
REST  . 

STACK:  SIGNAL  (NOW  SHORTER),  PREWHITENING  FILTER 

6  SETW  OFF  . 

SET  ANALYSIS  WINDOW  OFF  THE  SIGNAL  TRACE  . 

STACK:  SIGNAL,  PREWHITENING  FILTER 

7  WINDOW  COSINE  . 

WINDOW  THE  TOP  SIGNAL  WITH  A  COSINE  WINDOW  . 

STACK:  SIGNAL,  PREWHITENING  FILTER 

8  PAD  TO  1024  . 

INCREASE  THE  LENGTH  OF  THE  SIGNAL  TO  1024  POINTS  BY  ADDING  ZEROS  . 
NECESSARY  TO  MAKE  IT  COMPARABLE  TO  THE  PREWHITENING  FILTER. 

STACK:  SIGNAL,  PREWHITENING  FILTER 

9  FFT  . 

TAKE  THE  FFT  OF  THE  TOP  SIGNAL  AND  STORE  THE  REAL  AND  IMAGINARY 
PART  ON  TOP  OF  THE  STACK. 

STACK:  SPECTRUM  (REAL),  SPECTRUM  (IMAGINARY),  SIGNAL,  PREWHITENING 
FILTER 

10  SQR  2  . 

SQARE  THE  TWO  TOP  SPECTRA  . 

STACK:  SPECTRUM  (REAL),  SPECTRUM  (IMAGINARY),  SIGNAL,  P RE WHITENING 
FILTER 

11  ADD  . 

ADD  THE  TWO  TOP  SPECTRA  . 

STACK:  ENERGY  DENSITY  SPECTRUM,  SIGNAL,  PREWHITENING  FILTER 


Fig.  VII. 5. 3  PW  spectra  (page  2  of  3) 


12  SWAP  . 

SWAP  THE  TWO  FIRST  SIGNALS  . 

STACK:  SIGNAL,  ENERGY  DENSITY  SPECTRUM  ,  PREWHITENING  FILTER 

13  FLUSH  . 

FLUSH  OFF  THE  TOP  SIGNAL  . 

STACK:  ENERGY  DENSITY  SPECTRUM  ,  PREWHITENING  FILTER 

14  SWAP  . 

SWAP  THE  TWO  FIRST  SIGNALS  . 

STACK:  PREWHITENING  FILTER,  ENERGY  DENSITY  SPECTRUM 

15  FFT  . 

TAKE  THE  FFT  OF  THE  TOP  SIGNAL  AND  STOKE  THE  REAL  AND  IMAGINARY 
PART  ON  TOP  OF  THE  STACK. 

STACK:  SPECTRUM  (REAL),  SPECTRUM  (IMAGINARY),  PREWHITENING  FILTER, 
ENERGY  DENSITY  SPECTRUM 

16  SQR  2  . 

SQARE  THE  TWO  TOP  SPECTRA  . 

STACK:  SPECTRUM  (REAL),  SPECTRUM  (IMAGINARY),  PREWHITENING  FILTER, 
ENERGY  DENSITY  SPECTRUM 


17  ADD  . 

ADD  THE  TWO  TOP  SPECTRA  . 

STACK:  ENERGY  DENSITY  SPECTRUM  (PREWHITENING  FILTER),  PREWHITENING 
FILTER,  ENERGY  DENSITY  SPECTRUM 


18  SWAP  . 

SWAP  THE  TWO  FIRST  SIGNALS  . 

STACK:  PREWHITENING  FILTER,  ENERGY  DENSITY  SPECTRUM  (PREWHITENING 
FILTER),  ENERGY  DENSITY  SPECTRUM 


19  FLUSH  . 

FLUSH  OFF  THE  TOP  SIGNAL  . 

STACK:  ENERGY  DENSITY  SPECTRUM  (PREWHITENING  FILTER),  ENERGY  DENSITY 
SPECTRUM 


20  SWAP  . 

SWAP  THE  TWO  FIRST  SIGNALS  . 

STACK:  ENERGY  DENSITY  SPECTRUM,  ENERGY  DENSITY  SPECTRUM 
(PREWHITENING  FILTER) 


21  DIVIDE  . 

DIVIDE  THE  TWO  TOP  SPECTRA  POINT  BY  POINT  TO  ACCOUNT  FOR  THE  EFFECT 
OF  PREWHITENING. 

STACK:  ENERGY  DENSITY  SPECTRUM 

22  TERMINAL  . 

RETURN  FROM  THE  MACRO  MODE  . 
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VII. 6 _ Analysis  of  teleseismic  P  coda  using  NORSAR  and  NORESS 

The  teleseismic  P  coda  consists  of  the  presumably  scattered  energy 
that  arrives  in  the  first  few  minutes  following  first  P  from  tele¬ 
seismic  events.  It  has  been  suggested  that  it  could  be  used  as  a  yield 
estimator  (see,  for  example,  Ringdal,  1983),  and  work  reported  here 
also  suggests  that  it  may  also  be  useful  as  a  discriminant  for  crustal 
events.  This  study  examines  the  problem  of  the  teleseismic  P  coda  from 
the  standpoint  of  determining  the  place  where  the  waves  in  the  coda 
are  scattered,  with  the  aim  of  producing  a  quantitative  model  of  the 
coda. 

The  "starting  model"  is  shown  in  Fig.  VII. 6.1.  The  P  coda  is  con¬ 
sidered  to  be  produced  by  scattering  of  Lg  to  P  in  the  crust  near  the 
source  and  by  scattering  of  P  to  Lg  near  the  receiver.  The  power 
spectrum  Pc  of  the  coda  as  a  function  of  time  t  then  has  the  form,  for 
single  scattering, 

Pc(f,t)  -  |A|2  •  G  •  exp[-2xft/Q]  (1) 

where  f  is  frequency,  A(f)  is  the  Fourier  transform  of  the  first  P 
pulse,  G  is  a  scattering  parameter  called  the  turbidity  and  Q  is  the 
Q.  The  parameters  G  and  Q  will  refer  only  to  the  receiver  site  for  the 
case  of  a  deep-focus  event  and  will  depend  on  both  the  source  and 
receiver  sites  for  the  case  of  explosions  and  shallow  earthquakes. 

Fig.  VII. 6. 2  shows  a  plot  of  Fourier  amplitude  against  time  for  a 
Semlpalatinsk  explosion  for  a  frequency  of  3.6  Hz.  At  this  high  fre¬ 
quency  the  exponential  decay  model  of  Eq.  (1)  seems  to  fit  the  data 
well  for  the  time  range  20  sec  to  200  sec  after  first  P.  For  lower 
frequencies,  the  influence  of  PP  is  important  (King  et  al,  1975). 
Values  for  Q  of  850  ±  50  (deep-focus  events),  1200  ±  150  (Serai- 
palatinsk  explosions)  and  1200  ±  250  (Semlpalatinsk  crustal  earth¬ 
quake,  March  20,  1976;  Pooley  et  al,  1983)  have  been  obtained.  These 
values  suggest  Q  at  both  NORSAR  and  Semlpalatinsk  is  high.  Values  of  G 
of  about  10”^  km”*  are  found  for  both  deep-focus  events  and  Semi- 


palatinsk  explosions,  suggesting  that  the  size  of  the  coda  generated 
near  Semipalatinsk  is  at  most  of  the  same  order  as  that  generated  near 
NORSAR.  On  the  other  hand,  the  crustal  earthquake  near  the  Semi¬ 
palatinsk  test  site  gives  a  value  for  G  of  around  5  x  10“^  km-*,  con¬ 
siderably  larger  than  the  explosion  value.  This  can  be  explained  by 
efficient  generation  of  Lg  by  a  dislocation  source  followed  by  scat¬ 
tering  of  Lg  to  teleseismic  P. 

To  try  and  determine  the  relative  contribution  of  scattering  near  the 
source,  near  the  receiver  and  from  any  other  mechanism,  analysis  of 
NORESS  data  has  begun.  NORESS  is  uniquely  capable  of  determining 
crustal  phase  velocities  (3.5  -  8  km/sec)  and  azimuths  at  high 
frequencies.  From  Fig.  V1I.6.1,  the  model  predicts  that  a  component  of 
the  coda  should  have  a  phase  velocity  and  azimuth  appropriate  to  tele- 
seismic  travel  from  the  source,  and  another  component  having  phase 
velocities  appropriate  to  Lg  (3.5  -  4.5  km/sec)  and  random  azimuth. 
Fig.  VII. 6. 3  is  a  contour  plot  of  a  wavenumber  spectrum  of  a  5-second 
window  in  the  coda  of  a  Semipalatinsk  explosion,  5  seconds  after  P 
onset.  There  is  indeed  evidence  of  both  a  high  phase  velocity  peak  at 
the  appropriate  azimuth,  and  energy  at  Lg  velocities  off  azimuth. 

A.  Dainty,  Georgia  Inst, 
of  Technology 
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Fig.  VII. 6.1  Model  of  coda  (Top).  Generation  of  coda  by  Lg  +  P  scat¬ 
tering  near  source!  P  -*■  Lg  scattering  near  receiver 
(Middle).  Detail  of  Lg  >  P  scattering  near  source 
(Bottom).  Formulas  for  coda  power  spectrum  as  a  function 
of  time  in  time  range  20  -  200  sec  after  first  P.  See 
text  for  discussion. 


TIME  MINUTES 
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Fig.  VII. 6. 2  Fourier  amplitude  as  a  function  of  time  for  a  Semi- 

palatinsk  explosion,  April  22,  1978.  Line  indicating  fit 
of  equation  (1),  ticks  indicate  range  over  which  fit  is 
taken.  Analysis  frequency,  3.6  Hz.  Note  contamination  of 
coda  by  PP. 
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Fig.  VII. 6. 3  Contour  plot  of  wavenumber  spectrum,  KX  linear  wave- 

number  East,  KY  linear  wavenumber  North,  for  a  5-second 
coda  window  60  sec  after  1st  P.  Analysis  frequency  is 
3.6  Hz,  linear  wavenumber  is  (1/wavelength).  Circles  are 
at  phase  velocities  of  8  km/sec  (inner)  and  3.6  km/sec 
(outer). 


VII. 7  3-component  seismogram  analysis 

A  properly  equipped  seismograph  station  or  array  always  includes  3- 
component  instrumentation  for  the  very  simple  reason  that  the  seismic 
wavefield  comprises  vertical  and  horizontal  ground  motions  and  com¬ 
binations  thereof.  Seismologists  have  for  many  years  successfully 
exploited  the  information  potential  of  3-component  records  for  wave 
propagation  modelling,  retrieval  of  structural  information 
(tomography)  and  source  parameters,  but  these  efforts  have  mainly  been 
tied  to  the  low  frequency  part  of  the  wavefield.  Likewise,  many 
efforts  have  been  invested  in  extracting  similar  information  in  the 
high  frequency  range,  say  of  1-10  Hz,  but  in  this  case  less  success¬ 
fully  as  judged  from  current  literature.  The  reason  for  this  appears 
to  be  twofold:  1)  high  frequency  records  are  rather  complex  due  to 
scattering,  mode  conversions,  multipathing  and  similar  wave  propaga¬ 
tion  effects  and  ii)  the  analysis  techniques  used  fail  to  produce 
extracted  wavefield  parameters  in  an  easily  interpretable  format.  For 
example,  a  common  procedure  is  to  produce  particle  motion  plots 
reflecting  the  structure  in  the  wavefield,  but  since  such  plots 
generally  are  messy,  the  analysis  is  often  left  at  that. 

The  problem  addressed  in  this  note  is  that  of  a  new  approach  to 
extracting  parameter  characteristics,  the  wave  field  structure  on  the 
basis  of  a  priori  models  for  P,  S,  Love  and  Rayleigh  wave  particle 
motions.  Special  attention  has  been  given  to  the  problem  of  presenting 
results  in  an  easily  interpretable  manner  or  extracting  signal  parame¬ 
ters  convenient  for  a  wide  variety  of  research  applications.  These 
points  will  be  amply  demonstrated  in  the  Results  section. 

3-component  analyzing  technique  -  wavefield  modelling 

Any  approach  to  extracting  signal  parameters  from  a  seismogram  is 
based  on  certain  assumptions  or  models  regarding  the  signals  or  wave¬ 
lets  constituting  the  records  and  the  background  noise.  In  our  case 
these  are:  i)  noise  is  orthogonal  between  the  3-coraponents  (vertical, 
radial,  transverse;  or  Z,  E-W,  N-S)  and  ii)  particle  motion  is  as  derived 
from  classical  wave  theory.  The  analyzing  technique,  as  used  in  prac- 


Clce,  is  tied  to  the  theoretically  expected  particle  motion  covariance 
matrix  which  is  compared  to  that  observed  in  given  time  and  azimuth 
windows.  Then  using  a  max.  likelihood  formulation  the  primary  result 
is  simply  the  likelihood  that  a  wavelet  is  either  P,  S,  Love  or 
Rayleigh  against  the  null  hypothesis  of  being  noise.  Concurrent  with 
the  likelihood  estimation  is  that  of  extracting  the  axis  of  the  par¬ 
ticle  motion  ellipse  (linear  for  body  waves)  which  in  turn  can  be  con¬ 
verted  into  angle  of  approach  of  the  incoming  wavefront,  that  is, 
angle  of  incidence  (azimuth  is  already  known).  Details  on  this  par¬ 
ticular  particle  motion  analyzing  technique  is  presented  by  Christof- 
fersson  et  al,  1985.  To  summarize,  wavefield  information  extraction 
from  our  3-component  signal  analysis  procedure  is  as  follows: 

Likelihood  of  presence  (0.4  triggering  level)  of  P,  S,  Love 
and  Rayleigh  type  of  wavelets  as  a  function  of  azimtuh  (dAzi 
~  1  to  5  deg),  and  time  (dT  ~  1  sec,  updating  interval 
0.5  sec). 

Angle  of  incidence  for  triggered  body  wave  type  of  wavelets. 

As  mentioned,  special  attention  was  given  to  results  or  extracted 
wavefield  parameter  presentations  which  presently  are  in  the  following 
forms: 

i)  Likelihood  contouring  as  a  function  of  azimuth  (ranges  typi¬ 

cally  60-180  degrees)  and  time  (ranges  typically  20-50  sec). 

il)  Particle  motion  filtering  on  the  basis  of  the  estimated  like¬ 

lihood  function,  that  is,  simply  weighting  the  original 
records  with  the  likelihood  function  and  confining  to  an  azi¬ 
muth  window  of  ±  30  to  ±  40  deg  relative  "true"  azimuth. 


iii)  Angle  of  Incidence  estimates  converted  to  apparent  phase 
velocity  (so  far  for  P  only)  via  Herrin  tables,  and  then 
plotted  as  a  function  of  time  and  azimuth. 

lv)  Azimuth  and  incident  angle  of  the  very  first  P  wavelet  trans¬ 

formed  to  estimates  of  epicenter  coordinates. 

Ue  remark  that  presently  our  work  on  3-component  seismogram  analysis 
is  just  in  a  rather  preliminary  stage;  the  use  of  other  types  of  Mis¬ 
estimates  has  not  been  explored,  nor  the  use  of  more  complex  particle 
motion  models.  Other  types  of  problems  are  those  of  better  handling 
interference  phenomena  typical  of  the  Lg  wavetraln  and  S-wave 
splitting,  but  a  requirement  here  seems  to  be  access  to  S-component 
instrumentation  (1  vertical  and  4  horizontal).  On  the  application 
side,  we  have  not  explored  the  possibility  of  coda  decomposing  in 
deterministic  and  random  scattering  contributions,  and  perhaps  most 
challenging  that  of  what  we  term  earth  fingerprinting.  With  this  is 
meant  that  the  large  number  of  signal  parameters  extracted  may  exhibit 
stationary  patterns  reflecting  structural  heterogeneities  charac¬ 
teristic  of  specific  site  and  source  regions. 

It  may  also  be  appropriate  with  a  few  words  on  how  our  technique  com¬ 
pares  with  others  tied  to  analysis  of  3-component  data.  The  feature  in 
common  is  that  of  estimating  the  axis  of  the  particle  motion  ellipse, 
but  to  our  knowledge  nobody  has  tried  to  include  model  fit  in  probabi¬ 
listic  terms  nor  to  use  sliding  time/azimuth  windows  permitting  a 
rather  comprehensive  decomposition  of  the  whole  recorded  wavetraln. 

The  preliminary  results  to  be  presented  stem  from  analysis  of  the 
3-component  recordings  from  the  new  NOKESS  array,  and  obviously  some 
comparison  has  to  be  made  between  outcome  of  our  decomposition  tech¬ 
nique  and  similar  results  obtained  by  f-k  analysis  using  all  23  ver¬ 
tical  components  of  the  array.  Without  going  into  detail  here,  it 
suffices  to  state  that  3-component  results  (single  site)  compare 
favorably  to  the  f-k  results;  in  fact  quite  often  they  do  better  in 
terms  of  improved  time  and  azimuth  resolutions.  The  reason  for  this  is 
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that  f-k  analysis  requires  longer  time  windows  (a  minimum  of  2.5  sec 
is  used)  to  account  for  move-outs  across  the  array,  and  in  case  of 
interfering  signals  "averaged"  results  are  produced.  Not  to  forget, 
the  consistency  and  stationarity  of  particle  ground  motion  even  for 
short  wavelets  are  at  least  for  us  unexpectedly  good. 

Preliminary  results 

It  has  taken  considerable  time  to  program  our  3-component  analyzing 
technique;  the  first  rough  version  became  operational  quite  recently, 
so  comprehensive  analysis  of  many  event  recordings  has  not  been 
completed  yet.  However,  preliminary  results  from  a  few  rather  typical 
event  recordings  mostly  from  the  NOkESS  3-component  station  C2  will  be 
presented  in  the  following,  so  as  to  illuminate  in  our  opinion  the 
strength  and  potential  of  3-component  analyzing  techniques.  Results 
are  presented  event-wise  (in  appendices)  using  unfiltered  data  as 
standard  recursive  Butterworth  filtering  produces  severe  phase 
distortions.  It  should  be  added  that  so  far  mostly  P-modelling  has 
been  attempted. 

Case  1,  Event  No.  8A301:  Caspian  Sea/W.  Kazakh.  Feature:  travel  time 
triplication(s) . 

Case  2,  Event  No.  84327:  Local  explosion.  Feature:  complex  local 
record;  energy  migration. 

Case  3,  Event  No.  84328:  E.  Kazakh  event.  Feature:  triplication 
of  wave  train  ? 


Case  4,  Event  No.  85041:  Semipalatinsk  event.  Feature:  apparent 
velocity  analysis  of  the  first  25  sec  of  the  record. 


Case  5 ,  Event  No.  85060:  "Leningrad"  event.  Feature:  this  is 
actually  a  double  event  easily  separated  by  3-component 
analysis. 

Case  6,  Event  No.  85081:  Finland/USSR  event.  Feature:  worst  case 
event  with  exceptionally  poor  SNR,  but  3-component  analysis 
still  seems  to  work. 

Concluding  remarks 

The  preliminary  results  obtained  from  our  variant  of  particle  motion 
analysis  of  high-quality  3-component  records  from  the  new  NORESS  array 
in  Norway  are  considered  very  encouraging  and  thus  justify  research  on 
a  broader  scale  into  this  kind  of  problems.  Methodological  problems 
are  those  of  more  eficient  ML  estimators,  proper  parameter  settings 
including  shorter  updating  rates,  and  potential  for  real-time  event 
detection.  In  the  latter  case,  fast  algorithms  are  clearly  needed. 

In  fact,  a  solution  of  the  latter  problems  has  now  been  obtained,  that 
is,  the  probability  of  P  wave  presence  in  the  records  estimated  inde¬ 
pendently  of  azimuth.  However,  the  associated  calculations  are  rather 
demanding  and  Indeed  not  considered  justified.  Our  recommendation  is 
to  use  a  Walsh  detector  but  using  a  low  threshold  (high  false  alarm 
rate)  and  then  introducing  3-comp,  analysis  for  testing  precursor  of  P 
waves  in  an  off-line  mode. 

Surprisingly,  epicenter  location  capabilities  appear  to  be  relatively 
good,  azimuth  estimates  so  far  seldom  exceeded  ±  5  deg  (often  around 
±  2  deg),  while  distance  is  more  problematic  unless  secondary  phases 
(triplications)  are  clearly  identified.  In  the  same  way  as  the  azimuth 
angle  is  measured  in  the  3-comp,  analysis,  we  can  measure  similarly 
for  the  angle  of  incidence  in  the  vertical  plane  and  then  convert  to 
epicentral  distance  via  standard  travel  time  tables.  Experiments  with 
a  Semipalatinsk  event  gave  worst-case  distance  error  of  5  deg  (see 
also  Sec.  V11.8).  Another  interesting  aspect  of  3-comp,  seismogram 
analysis  is  that  of  identifying  pP  and  thus  significantly  improving 


focal  depth  estimates  which  besides  are  highly  diagnostic  of  source 
type. 
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CASE  1 


Event  No.  84301: 


Caspian  Sea/U.  Kazakh 
PDE  azi:  ~  106°,  A  ~  25° 

NORESS  azi:  -  110°,  vel.  -  12.4  km/s 


Note:  Triplication  features 


Fig.  A:  Original  unfiltered  data 
Fig.  B:  Probability  of  P-wave  presence 
Fig.  C:  Probability  function  of  Fig.  B  used  as  a  filter 
Fig.  D:  King  &  Calcagnile  travel  time  model;  epicenter  distance 
marked 


Comments:  The  particle  motion  filter  "produces”  phase  arrivals  in 
very  good  agreement  with  expected  travel  time  curves.  The  surplus 
arrival  may  be  associated  with  a  hypothesized  discontinuity  at  ca 
550  km  depth.  Epicenter  locations  better  than  1  deg  are  deemed 
feasible  in  such  cases,  as  triplications  when  identified  give  very 
precise  distance  estimates. 
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Event  No.  84327: 


Local  explosion. 
NORESS:  Azi  -  "south" 


Fig.  A:  Original  data 

Fig.  B:  Probability  of  P-wave  presence;  note  trend  of  migrations  in 
pattern 

Fig.  C:  Particle  motion  filtered  output 


Comments:  At  a  later  stage  filters  based  on  S  (SV  &  SH),  SV,  Love  and 
Rayleigh  presence  will  be  introduced.  The  off-azimuth  triggering  at 
about  260-300  deg  may  reflect  sort  of  side  lobe  effects.  Anyway, 
3-component  filtering  of  complex  local  records  seems  clearly  to  be 
feasible. 
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CASE  3 


Event  No.  84328: 


Eastern  Kazakh 


NORESS  gives  location  azi  *  79°,  vel.  =  13.8  km/s 


Fig.  A:  Original,  unfiltered  data 

Fig.  B:  Probability  of  P-wave  presence;  note  triggering  before  main 
arrival 

Fig.  C:  Filtered  data;  indications  of  triplication  plus  an  unidenti¬ 
fied  precursory  arrival 


Comment:  With  practical  experience.  Fig.  C  filtered  output  may  be 
associated  with  a  specific  siting  area  and/or  distance  range.  Other 
problem,  false  alarm  rate  of  probability  triggering. 
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CASE  4 


Event  No.  85041 


Semipalatinsk 


PDE  azi  -  74°,  A  ~  38° 


NORESS  azi  -  79°,  vel.  13.8  km/s 


Note:  Exceptionally  clear  records 


Fig.  A:  Original,  unfiltered  records 

Fig.  B:  Probability  of  P-wave  presence;  parameter  setting  extremely 
restricted  so  as  to  exclude  secondary  arrivals 

Fig.  C:  Particle  motion  filtered  output 

Fig.  0:  Apparent  velocity  of  detected  P-wavelets  under  less  severe 
triggering  conditions  -  time  interval  25  sec  after  first 
onsets 


Comments:  Fig.  D  demonstrates  feasiblity  of  signal  coda  decom¬ 
position;  results  in  good  agreement  with  f-k  analysis  using  all  NORESS 
vertical  sensors  (Or.  A.  Dainty,  personal  communication). 
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CASE  5 


Event  No.  85060 


Leningrad  event;  azi.  ~  87°,  A  ~  8.8° 


NORESS:  azi.  ~  87°,  vel.  -  12.1  km/a 


Note:  This  is  actually  a  double  event,  the  first  arriving  at  approx. 
120  sec,  azi.  ~  200°,  and  "Leningrad"  at  127  sec,  azi.  ~  90°. 
Both  seen  in  the  3-component  analyzing  results. 


Fig.  A:  Original,  unfiltered  traces;  the  two  arrival  marked  1  and 
2,  respectively 

Fig.  B:  Probability  of  P-wave  presence 
Fig.  C:  Particle  motion  filtered  output 


Comments:  Our  3-component  analyzing  technique  appears  to  have  some 
potential  for  real-time  event  detection. 


CASE  6 


Event  No.  85081 

Finnish  epicenter  location  at  Finland/USSR  Border  gives  azi.  ~  71°, 
A  ~  9.8° 

NORESS:  azi.  ~  90°,  vel.  -  16.7  km/s 


Note:  The  unfiltered  data  used  in  our  analysis  have  exceptionally 

poor  SNR;  phase  arrival (s)  relatively  clearly  visible  in  full 
NORESS  record  display. 


Fig.  A:  Original,  unfiltered  data 
Fig.  B:  Probability  of  P-wave  presence 
Fig.  C:  Particle  motion  filtered  record 


Comments:  P-wave  particle  motion  structure  appears  to  be  preserved  to 
some  extent  even  for  poor  SNRs. 
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VII. 8  Epicenter  location  scenarios 

Traditionally,  P-wave  arrival  time  observations  from  4  stations  are 
required  for  the  estimation  of  the  four  focal  parameters  epicenter 
location  (o,  X),  focal  depth  and  event  origin  time.  The  computational 
techniques  in  use  are  well  documented  in  the  scientific  literature, 
although  an  occasional  severe  practical  problem  is  that  of  proper 
phase  association.  The  latter  difficulty  stems  from  inadequate 
grouping  of  event  P  reportings  when  several  events  "overlap"  in  a 
given  time  window.  Such  ambiguities  are  easily  resolved  when 
array/3-comp.  station  estimates  are  available.  In  this  section  three 
scenarios  for  event  location  will  be  discussed,  namely:  i)  single  sta¬ 
tion  (3-comp.)  array  event  location;  ii)  two  station/array  event  loca¬ 
tion;  and  iii)  precise  event  location  using  a  regional  network  of  high 
quality  stations. 

i)  Single  station/array  event  location 

The  problem  is  illustrated  In  Fig.  VII. 8. la  and  VII. 8. 2,  and  obviously 
the  accuracy  in  the  epicenter  location  estimate  reflects  precision  in 
azimuth  and  distance  estimates.  As  demonstrated  in  Section  VII. 7,  a  single 
3-comp,  station  can  provide  azimuth  estimates  within  5  degrees  or  even 
within  2  degrees  in  good  cases.  Distance  estimates  are  more  problema¬ 
tic  unless  secondary  phases  confidently  can  be  identified.  This  is 
indeed  feasible  at  local  (Pn,  Sn,  Lg,  etc.)  and  regional  (400  km  and 
650  km  travel  time  triplication  curves)  distances,  while  slowness 
estimates  are  not  too  accurate  at  teleseismic  distances.  Using  3-comp, 
recordings,  the  measured  angle  of  incidence  can  provide  distance  esti¬ 
mates  which  appear  to  be  more  precise  than  those  stemming  from 
slowness  measurements  say  via  F-k  analysis  of  dat  afrom  small  aperture 
arrays.  Another  advantage  with  3-corap.  data  analysis  is  that  phase 
identification  becomes  more  objective,  and  may  even  provide  a  means 
for  focal  depth  estimation  via  "detection"  of  the  pP-phase. 


i 

I 
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ii)  Two  station/array  event  location 

The  problem  is  illustrated  in  Fig.  VII. 8. lb  and  VII. 8. 3,  and  the  most 
important  parameters  needed  are  azimuths  and  relative  P  arrival  times 
as  a  substitute  epicentral  distance  estimates.  With  reasonable 
geometries,  epicenter  location  should  be  within  a  few  hundred  kilo¬ 
meters  or  even  less. 

Comments.  Gj0ystdal  et  al  (1973)  have  extensively  dealt  with  one-  and 
two-array  event  location  capabilities  and  theic  techniques  naturally 
apply  to  azimuth  and  distance  estimates  extracted  from  the  3-corap. 
station  data.  Naturally,  the  accuracy  in  the  epicenter  locations  made 
depends  on  geometry  (see  Fig.  VI1.8.1),  as  wavefront  angles  scale  dif¬ 
ferently  with  distances.  An  interesting  problem  here,  now  under 
investigation,  is  whether  the  pattern  of  secondary  arriving  phases  as 
derived  via  our  new  3-comp,  data  analysis  can  be  used  for  refining 
epicenter  distance  estimates. 

iii)  Use  of  regional  networks  for  event  locations 

In  the  context  of  a  potential  comprehensive  test  ban  treaty,  precise 
event  locations  are  a  necessity,  and  this  can  only  be  achieved  by 
having  access  to  recordings  from  a  network  of  regional  stations  (for 
small  events).  A  number  of  computer  schemes,  such  as  HYPO  75,  exist 
for  solving  this  kind  of  problems,  but  a  general  drawback  here  is  that 
the  crustal  model  used  is  considered  a  bit  crude.  For  example,  what  is 
the  velocity  distribution  in  the  crust?  Is  it  constant,  depth  depen¬ 
dent  only;  is  lateral  variation  significant  or  do  we  have  to  take 
potential  anisotropy  into  account?  A  related  problem  is  which  phases 
are  actually  observable  in  the  seismogram  and  what  are  their  respec¬ 
tive  ray  paths.  Do  we  actually  see  waves  refracted  from  the  Conrad 
discontinuity;  the  Pn  phase  may  be  a  head  wave  or  perhaps  better 
modelled  as  a  diving  wave?  Likewise,  reflections  from  Conrad  and/or 
Moho  if  observable  would  constrain  focal  depth  estimates.  Clearly, 
precise  event  locations  using  regional  recordings  are  much  dependent 
on  the  velocity  distribution  in  the  crust  and  lithosphere  and  our  abi¬ 
lity  to  model  it  properly,  and  equally  important,  extract  relevant 
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phase  information  from  the  seismogram.  Extensive  use  of  S-phases  is 
clearly  an  advantage  here  in  view  of  their  short  wave  lengths. 


A  long-term  research  project  is  now  under  way  at  NORSAR  for  solving 
these  kinds  of  problems;  the  principal  goals  are  those  of  Inverting 
available  travel  time  information  for  a  multitude  of  phases  in  order 
to  produce  a  tomographic  image  of  the  crust  and  lithosphere  structures 
and,  important,  jointly  extract  relevant  focal  parameters.  The 
research  strategy  here  is  outlined  in  some  more  detail  by  Gubbins 
(1985),  while  Christoff ersson  et  al  (1985)  have  presented  a  novel 
approach  to  both  quantitative  and  qualitative  extraction  of  signifi¬ 
cantly  more  phase  information  from  3-comp,  seismogram  recordings 
(e.g.,  see  Section  VII. 7). 

A.  Christoffersson, 

Uppsala  Univ. 

D.  Gubbins,  Cambridge 
Univ. 

E. S.  Husebye 
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VII. 8. 2  The  95  per  cent  confidence  ellipses  of  simulated  one- 

array  (or  one  3-corap.  station)  location  capabilities  for 
epicentral  distances  of  30°  and  80°.  Standard  deviations 
of  azimuth  and  slowness  are  1.5°  and  0.1U  sec/deg, 
respectively,  and  the  number  of  trials  is  200.  Note  how 
the  axis  of  the  error  ellipse  changes  with  distance 
(figure  after  Gj0ystdal  et  al,  1973). 
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Fig.  VII. 8. 3  Simulated  2-array  (or  two  3-comp,  stations)  event  loca¬ 
tions  at  45 . 3°N;  149. 3°E  for  NORSAR  and  L.\SA.  In  case  a) 
only  azimuth  observations  are  used  (o  *  2  deg),  while  in 
case  b)  further  constraints  are  added  by  using  P  arrival 
time  differences  (o  =*  2.0  sec). 


VII. 9  The  RST/RSTE  research  project  -  current  status 

The  RST/RSTE  (Remote  Seismic  Terminal  Enhanced)  project  is  aimed 
at  developing  an  intelligent  seismic  field  terminal  with  capabilities 
of  detecting  and  analyzing  events  in  an  automated  mode.  The  opera¬ 
tional  capabilities  would  be  roughly  equivalent  to  those  of  the  dedi¬ 
cated  data  center  for  a  large-aperture  array,  although  the  number  of 
channels  to  be  handled  would  be  far  less.  This  entails  a  multitude  of 
tasks  like  real-time  event  processing,  off-line  processing  of  detected 
events,  and  external  transfer  to  an  RST  of  raw  and  processed  signal 
data.  The  viability  and  realizability  of  the  RST/RSTE  concept  was 
demonstrated  in  Geneva  in  1982  (e.g.,  Husebye  and  Thoresen,  1984);  for 
a  schematic  view,  see  Fig.  VI I. 9.1. 


The  mentioned  operational  task,  say  for  8  short  period  seis¬ 
mometers,  can  be  "accomodated"  with  a  hardware  framework  of  mini  -or 
personal  computers  like  the  IBM  PC  AT.  The  corresponding  software 
implementation  requires  a  highly  efficient  multitasking  operating 
system,  as  well  as  a  high  level  programming  language.  Our  choice  was 
polyFORTH  II,  which  is  an  integrated  operating  system/high  level 
language  (FORTH).  The  detection  process  is  a  voting  one  tied  to 
STA/LTA  thresholding  on  individual  channels  -  a  Ualsh  detector  re¬ 
placement  is  under  consideration.  In  the  event  processor  detected 
signals  would  be  subjected  to  f-k  analysis,  and  a  suitable  phase 
picker  for  signal  onset  times.  Important,  in  view  of  the  very  recent 
3-comp,  seismogram  analysis  results  reported  in  Section  VII. 7,  a 
viable  alternative  here  would  be  to  replace  the  mini-array  with 
3-comp,  high-quality  instrumentation.  It  should  be  added  that  recent 
progress  in  fault  mechanism  studies  entail  that  reasonable  moment 
tesnor  solutions  are  feasible  in  many  cases,  even  when  single  3-comp, 
station  records  are  used  (prof.  A.  Dziewonski,  personal 
communication) . 


Seismological  motivations  for  RST/RSTE 

Theoretical  advances  in  seismology  cannot  be  sustained  unless 
researchers  have  easy  access  to  an  abundance  of  high-quality  seismic 
data  in  digital  form.  The  NORSAR  array  and  the  SRO  and  ASRO  seismic 
observatories  are  typical  examples  of  realizations  of  these  data 
collection  goals.  On  the  other  hand,  most  seismometers  in  operation 
are  tied  to  analog  recording  systems  (photographic  paper  to  be  changed 
daily).  In  other  cases  the  seismometer  outputs  are  transmitted  con¬ 
tinuously  say  via  dedicated  telephone  lines  to  a  data  center,  although 
only  a  few  per  cent  of  the  data  flow  contains  significant  event 
information.  Also,  the  use  of  small  arrays/3-comp.  stations  will  very 
much  facilitate  signal  analysis,  that  is,  Level  1  data  extraction, 
both  qualitatively  and  quantitatively.  The  obvious  thing  to  do  here  is 
naturally  to  decide  at  the  seismometer  site  which  part  of  the  recor¬ 
dings  are  useful,  and  in  this  way  greatly  reduce  excessive 
transmission  costs.  Likewise,  extraction  of  Level  1  data  in  the  field 
would  permit  much  faster  and  more  advanced  event  analysis  with  a  mini¬ 
mum  of  analyst  interaction  and  delays.  It  should  be  noted  that  if  mass 
data  storage  in  the  field  is  found  desirable  by  any  user,  this  can  be 
achieved  by  installing  minitape  recorders.  For  example,  makes 
currently  available  can  store  about  70  Mbytes  of  data  (cost  ~  $3000), 
while  most  recent  models  can  store  ca  130  Mbytes  (at  the  same  pri¬ 
ce)  which  is  equivalent  to  the  standard  9T  6250  BPI  computer  tape. 

We  feel  confident  that  our  hardware/software  designs  of  the  RST  and 
RSTE  to  be  detailed  in  the  next  chapters  will  accomplish  these  goals. 
Finally,  let  us  add  that  in  the  case  of  broad-band  seismograph  sta¬ 
tions  the  data  recording  is  primarily  of  interest  for  research  and  in 
such  cases  the  RFT/RSTE  would  serve  as  a  data  logger  system. 

RSTE  and  RST  software  developments 

Basically,  our  intention  is  to  make  the  RSTE  a  smart  peripheral  device 
intevfaceable  to  a  variety  of  host  computers  or  RSTs.  Its  main  task  is 
time  series  analysis  in  the  field  and  data  transfer  on  the  basis  of 
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intelligent  handling  of  outputs  from  a  variety  of  geophysical  instru¬ 
mentation,  like  small  arrays  of  short  period  seismometers,  hydrophone 
arrays  and  3-component  broadband  seismometers.  Naturally,  we  do  not 
intend  to  incorporate  all  possible  functions  and  tasks  in  one  RSTE, 
instead  aim  at  designing  a  flexible  software  framework  from  which  spe¬ 
cial  purpose  RSTE8  can  be  realized  relatively  easily.  For  example, 
seismic  data  recording  and  subsequent  analysis  are  real-time  multi¬ 
tasking  processes  which  motivated  the  choice  of  polyFORTH  (Whitney  and 
Conrad,  1983)  as  programming  language.  This  is  an  integrated  operating 
system/high  level  language  with  minimal  overhead  for  multitasking. 
Likewise,  the  choice  of  the  KERMIT  protocol  for  two-way  RST-RSTE  com¬ 
munication  was  motivated  by  its  extensive  use  in  universities  and 
research  institutions  (e.g.,  see  DaCruz  and  Catchings,  1984a, b). 

The  RSTE  functional  structure  as  realized  in  the  software  is 
shown  in  Fig.  Vll.9.2.  All  boxes  in  the  figure  represent  independent 
and  concurrent  processes,  each  having  a  well-defined  function  and 
interface.  These  concurrent  processes  are  executed  and  synchronized  by 
the  Round  Robin  multitasker  of  polyFORTH  on  a  single  processor.  The 
arrows  visualize  the  communication  between  processes,  with  the  ver¬ 
tical  line  in  the  middle  representing  a  "software  bus".  All  processes 
connected  to  this  software  bus  may  "talk"  to  one  another.  The  actual 
Implementation  of  inter-process  communication  is  done  using 
"mailboxes".  Each  process  has  an  associated  mailbox  where  other  pro¬ 
cesses  may  put  messages.  Details  on  the  RSTE  software  processes 
shown  in  Fig.  VII. 9. 3,  most  of  them  now  completed,  are  given  below. 

KERMIT 

This  process  implements  the  communications  channel  to  the  outside 
world,  using  a  modified/extended  version  of  the  KERMIT  protocol.  The 
extensions  include  having  KERMIT  running  permanently  in  server-mode, 
login  procedure,  long  packets  (up  to  233  bytes)  and  transfer  of  8-bit 
data  without  quoting.  File  transfers  are  handled  in  cooperation  with 
the  File  System.  Commands  received  through  KERMIT  are  handed  over  to 
the  "Command  Interpreter"  for  parsing.  The  KERMIT  protocol  is  imple- 


mented  as  a  "finite  state  machine",  and  an  example  of  state  descrip¬ 
tion  of  the  protocol  is  illustrated  in  Fig.  VII. 9. 4.  It  is  a  pleasure 
to  note  that  a  seismic  data  logger  under  development  at  Harvard  Uni¬ 
versity  will  feature  a  KERMIT  protocol  very  similar  to  ours  (J. 

Steim,  personal  communication). 

MODEM 

This  is  the  interface  between  KERMIT  and  the  hardware  used  for  com¬ 
munications.  MODEM  accepts  commands  (DIAL,  DISCONNECT,  INITIALIZE, 

...)  and  data-packets  from  KERMIT,  while  sending  messages  (incoming 
call,  carrier  lost,  etc.)  and  data-packets  back  to  KERMIT. 

KA 

KA  is  short  for  "KERMIT  Allocator",  and  is  needed  because  the  com¬ 
munications  channel  cannot  be  shared.  Every  process  wishing  to  use 
KERMIT  must  first  allocate  it  for  exclusive  use,  and  this  is  done 
through  requests  to  KA. 

COMMAND  INTERPRETER 

This  process  is  responsible  for  parsing  and  execution  of  commands 
received  through  KERMIT. 

FILE  SYSTEM 

The  function  of  this  process  is  obvious:  it  offers  all  other  processes 
necessary  file-handling  functions.  Multiple  open  files  are  allowed. 

The  file  system  can  be  implemented  using  one  or  both  of  RAM  (memory) 
and  disk  storage.  The  actual  implementation  is  transparent  to  the  rest 
of  the  system.  Standard  DOS  naming  conventions  are  used  with  up  to  8 
characters  name  and  3  characters  extension. 

CONFIG 

Activated  by  a  command  from  the  COMMAND  INTERPRETER,  this  process 
enables  remote  configuration  of  the  entire  RSTE-system.  This  includes 
setting  filter  parameters,  sampling  frequency,  detection  sensitivity, 


type  of  post-processing  to  be  used  and  so  on.  The  required  con¬ 
figuration  is  specified  in  a  file  that  has  been  prepared  transparently 
and  transferred  in  advance  by  the  controlling  RST. 

OPSYS 

OPSYS  is  nothing  more  than  a  collection  of  utilities  and  system  func¬ 
tions  needed  by  other  processes. 

EVENT  FILE  DISPATCHER 

This  process  handles  dispatching  (and  possible  queuing)  of  event  files 
that  have  been  declared  "interesting"  by  the  EVENT  post-processor.  The 
EVENT  FILE  DISPATCHER  has  available  a  phone  directory,  which  is  used 
to  dial  an  RST  and  transfer  the  file. 

WATCHDOG 

The  WATCHDOG  is  responsible  for  monitoring  the  state  of  the  system, 
reporting  and  logging  minor  errors.  If  this  process  after  crash  does 
not  restart  and  external  timer,  that  timer  will  automatically  restart 
the  system. 

EVENT  POSTPROC 

The  EVENT  post-processor,  one  of  the  major  RSTE  modules,  consists  of 
various  analysis  programs  that  are  run  in  background  on  data  stored  by 
the  DATAFLOW  PROCESS.  These  are  data  blocks  which  are  part  of  the  out¬ 
put  from  the  real-time  event  detection  operation.  The  principal  task 
here  is  to  provide  a  listing  of  detected  events  in  terms  of  extracted 
Level  1  parameters  and  rough  epicenter  locations. 

TIME 

TIME  is  just  an  interfae  to  the  available  real-time  hardware,  offering 
precise  absolute  and  relative  time  to  all  processes. 


TEST  &  CALIBRATION 

This  process,  together  with  some  external  hardware  (oscillator),  pro¬ 
vides  the  ability  to  calibrate  seismometers. 

DATAFLOW  PROCESS 

This  is  the  critical  module  in  the  RSTE  handling  system.  The  seismometer 
data  stream  is  subdivided  into  many  small  concurrent  processes,  each  imple¬ 
menting  a  part  of  the  data-flow  like  the  multiplexed  digitizer,  data  win¬ 
dowing  and  filtering,  etc. 

Comments.  The  RSTE  software  system  is  now  nearing  completion  with  the 
EVENT  POSTPROC  being  a  major  exception  here.  The  reason  for  this  is  that  we 
are  currently  somewhat  ambivalent  about  whether  the  "best"  solution  would 
be  a  miniarray  or  a  high-quality  3-comp,  installation.  The  latter  alter¬ 
native  is  rather  attractive  not  only  as  regards  decomposition  of  complex 
recordings  but  also  for  simpler  field  logistics  (avoiding  stretching  of 
cables,  etc.). 

RST  functional  structure 

The  RST  functional  structure,  shown  in  Fig.  9,  is  grossly  similar  to 
that  of  the  RSTE.  Of  course,  the  RST  has  no  DATAFLOW  process  and  no 
EVENTFILE  DISPATCHER,  while  the  KERMIT/MODEM/KA  and  FILE  SYSTEM  per¬ 
form  the  same  functions  as  in  the  RSTE. 

There  is  little  concurrency  in  the  RST.  Only  two  tasks  are  sharing  the 
processor:  KERMIT/FILE  SYSTEM  and  the  COMMAND  INTERPRETER.  Boxes  to 
the  right  of  the  "software  bus"  denote  subprograms  (modules),  to  be 
outlined  below,  that  are  activated  one  at  a  time  under  control  of  the 
COMMAND  INTERPRETER,  which  is  responsible  for  the  user  interface. 

Signal  analysis  and  displays 

This  module  enables  the  user  to  do  different  kinds  of  analysis  on  data 
received  from  RSTEs.  This  includes: 


Display  of  traces  with  zoora/pan  on  both  color  graphics  screen  and 
printer/plotter.  Use  of  joy  stick/mouse  would  enable  an  analyst  to 
re-estimate  onset  time,  signal  amplitude  and  period,  etc. 

-  Spectral,  f-k  plots  and/or  3-comp,  wave  train  decomposition 

-  Joint  display  of  traces  from  more  than  one  RSTE. 

RSTE  configuration  utilities 

This  module  is  not  required  to  receive/analyze  event  files,  but 
instead  is  aimed  at  easy  preparation  of  configuration  files,  for 
use  in  the  remote  RSTE. 

Device-dr iver s 

Modules  named  "PRINTER",  "PEN  PLOTTER"  and  "GRAPHICS  DISPLAY"  are 
simply  device-drivers  used  by  the  various  display  functions  of 
ANALYSIS. 

Concluding  remarks 

Our  RST/RSTE  concepts  have  evolved  from  what  we  may  term  Personal 
Seismometry  Now,  as  demonstrated  In  Geneva  in  1982  (Husebye  and 
Thoresen,  1984)  to  "Seismic  Arrays  for  Everyone"  as  sketched  by  Huse¬ 
bye  et  al  (1985).  The  motivation  for  undertaking  this  project  is 
simply  the  obvious  need  for  the  selsmological  community  at  large  to 
have  easy  access  to  high-quality  short  period  seismic  data  and  without 
literally  being  ruined  in  the  process.  Similar  concerns  are  shared  as 
regards  broadband  seismometry  which  has  motivated  the  establishment  of 
regional  and  global  broadband  station  networks. 

The  RST/RSTE  system  is  modular  and  flexible,  and  this  can  be 
modified/extended  to  meet  special  user  modifications.  An  important 
step  in  this  direction  is  the  planned  replacement  of  the  present  RSTE 
microcomputer  OMNIBYTE  with  an  IBM  PC  AT.  In  the  extreme  this  would 


entail  RST/RSTE  to  serve  either  as  a  data  logger  or  at  the  other  end 
as  a  seismological  analyst  work  station. 


E.S.  Husebye 
E.  Thoresen,  Consultant 
J.A.  Mathiesen 
l.C.  Myrvold 
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Fig.  VII. 9. 2  The  RSTE  functional  structure.  Details  on  functional 
tasks  are  given  in  the  text. 


n  KERMIT  communication  protocol  state  transmission  diagrams  illustrating 
the  RSTE. 


VII. 10.1  Crustal  structure  in  the  NORSAR  array  siting  area  derived 
from  gravity  observations 

The  NORSAR  array  is  sited  at  the  outskirts  of  the  paleo  Oslo  Rift/ 

Oslo  Grabea  system,  which  is  geophysically  characterized  by  a  pro¬ 
nounced  gravity  high.  These  data  may  provide  details  on  the  crustal 
structure  in  the  NORSAR  area  as  the  array's  seismic  recordings, 
despite  much  research  ingenuity,  have  proved  to  be  essentially 
transparent  to  local  heterogeneities  due  to  near-vertical  angles  of 
incidence  of  the  teleseismic  waves  (e.g.,  see  Aki  et  al  1977; 

Christof fersson  and  Husebye,  1979;  Haddon  and  Husebye,  1979;  Troitskiy 
et  al,  1981;  and  Tompson  and  Gubblns,  1981).  Crustal  profiling  surveys 
and  local  earthquake  recordings  are  not  particularly  helpful  in  this 
respect  and  the  non-uniqueness  inherent  in  refraction  profiling  types 
of  surveys  (no  cross-sampling)  remain  a  problem  (e.g.,  see  Cassell  et 
al,  1983;  Mereu  et  al,  1983). 

Clearly,  crustal  structural  details  are  important  for  an  improved 
understanding  of  complexities  in  high  quality  seismic  records  like 
those  now  provided  by  the  new  NORESS  array,  and  as  a  first  step 
towards  solving  this  kind  of  wave  propagation/scattering  problems  we 
have  attempted  to  exploit  the  structural  information  implicitly  con¬ 
tained  in  the  easily  available  Oslo  Graben  gravity  data. 

Tectonic  setting  and  the  Oslo  Graben  residual  gravity  field 

The  tectonic  setting  of  the  general  Oslo  Graben  area  is  shown  in 
Fig.  VII. 10.1  and  clearly  implies  that  its  evolution  has  been  rather 
complex.  However,  on  the  other  hand  the  geophysical  imprints  of  past 
orogens,  etc.,  are  not  preserved  indefinitely,  and  today  the  Oslo  Gra¬ 
ben  gravity  high  appears  to  be  the  outstanding  feature  reflecting  past 
taphrogenesis.  The  gravity  data  coverage  of  the  general  graben  area  is 
relatively  dense,  altogether  more  than  5300  measurements  have  been 
carried  out  (e.g.,  see  Ramberg,  1976).  His  residual  gravity  field  is 


shown  in  Fig.  VII. 10. 2,  and  this  map  formed  the  basis  for  our  rein¬ 
terpretation  (inversion)  of  the  gravity  data  pertaining  to  the  Oslo 
Graben  area. 

For  the  sake  of  completeness  it  should  be  mentioned  that  Ramberg 
attributed  the  observed  anomalies  to  intrusion  of  relatively  dense 
material  at  the  base  of  the  crust  (Moho  upwarp)  while  Husebye  et  al 
(1978)  using  Parker's  (1975)  Ideal-Body  Concept  gave  a  mid-crust  loca¬ 
tion  of  the  intrusive  body.  Parker's  technique  is  very  different  from 
standard  approaches  as  in  order  to  reduce  the  non-uniqueness  inherent 
in  gravity  interpretations,  information  extraction  is  optimized  with 
respect  to  either  density  contrast  or  distance  to  causative  body. 

Linear  density  inversion 

The  gravity  field,  being  a  potential  field,  invites  ambiguity  in 
interpretations  which  was  duly  recognized  at  a  relatively  early  stage 
of  exploration.  However,  techniques  differ  in  their  ability  to 
reduce  the  number  of  plausible  solutions  without  Invoking  a  priori 
information  about  the  heterogeneous  body  such  as  density  contrasts  and 
geological  structural  boundaries  (outcrops).  Thus,  inverse  methods  are 
favored  because  they  require  a  minimum  number  of  subjective  judgements 
on  the  part  of  the  Interpreter. 

The  physical  model  used  (see  Fig.  VII. 10.3)  is  of  the  2$-D  kind  and 
consists  of  N  block-formed  prisms  of  200  km  extent  matching  the 
approximate  length  of  the  Oslo  Graben.  For  the  i-th  prism  the  gravita¬ 
tional  attraction  gjj  at  the  j-th  observation  point  is  a  linear  func¬ 
tion  of  density  p^,  that  is: 

.  -  1 _ 

8ij<Pi>  *  YPi  J  z  ’  ro  -  dV  -  aijPi 

V  I r-r0| 


(1) 
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where  y  “  universal  gravitational  constant,  and  the  volume  Integral 
equals  the  geometrical  weighting  funtion  ajj  relating  volume  of  the 
prism  and  its  distance  from  the  observation  point. 

The  theoretical  gravity  Gjfc  at  N  observation  points  is: 

>  M 

G/  -  Fj(p)  «^gij(Pi)  j  -  1,N  (2) 

To  estimate  a  density  distribution,  we  note  that  by  substituting  the 
"true"  density  values,  px>  into  eq.  (2),  we  get 


Gjt  -  Fj(pT)  -  Gj°  j  -  1,N  (3) 


where  Gj°  is  the  observed  gravity  at  point  j.  Since  Fj  is  linear  in 
density  and  noting  that  our  initial  density  model  p0  is  related  to  p? 
by 

>  ■*  + 

pT  -  P0  +  Ap  (4) 

we  achieve 

+  ■»  ♦ 

Fj(p0+Ap)  -  Fj(p0)  +  Fj(Ap)  -  Gj°  j  -  1,N  (5) 
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Eq.  (1),  (2)  and  (5)  combine  Co  yield 


Gj°  •  Fj(p0)  -  Fj(Ap) 


M 

I  *ijApi 
i-1 


j  "  1  »M 


Using  matrix  notation, 


(6) 


AG  -  G°  -  Gt  -  AAp  (7) 

A  least  squares  solution  to  eq.  (7)  for  Ap  is  obtained  by  solving  the 
normalized  system: 

ATA  Ap  “  A-^  AG  (8) 

where  T  denotes  transpose.  Now,  A  can  be  decomposed  as: 

A  -  U  A  VT  (9) 

Inserting  this  in  eq.  (8)  and  taking  the  generalized  inverse  of  ATA, 
the  estimate  of  the  correction  vector  becomes: 

Ap  -  V  A_1Ut  AG  -  H  AG  (10) 

The  Ap  estimation  procedure  is  strictly  linear  since  the  A-matrix 
remains  invariant. 

Important  by-products  of  the  inversion  procedure  are  the  resolution 
and  the  information  density  matrices,  defined  respectively  as: 

R  -  HA  -  V  VT 

S  -  AH  -  U  UT 


For  details  on  techniques  used  for  solving  the  eq.  (7)  type  of  linear 
equations,  reference  is  made  to  Backus  and  Gilbert  (1968).  We  want  to 
remark  that  the  presence  of  small  eigenvalues  may  cause  the  solution 
to  oscillate  markedly,  so  a  lower  cutoff  value  has  to  be  selected.  In 
practice,  this  was  achieved  by  retaining  the  variance  of  the  parameter 
estimate  relatively  small  without  sacrificing  too  much  of  the  resolu¬ 
tion  which  amounts  to  eliminating  the  influence  of  short  wavelength 
noise  in  the  data.  Added  flexibility  in  using  the  above  inversion 
scheme  is  often  desired  as  the  initial  systems  of  prisms  (or  blocks) 
may  at  a  later  stage  in  the  analysis  prove  to  be  somewhat  awkward.  For 
example,  it  may  be  desirable  to  reduce  the  number  of  unknowns  by 
merging  prisms,  or  to  build  in  a  priori  geological  information  about 
the  likely  shape  of  the  anomalous  body.  Christoff ersson  and  Husebye 
(1979)  have  demonstrated  a  simple  scheme  for  Incorporating  such  modi¬ 
fications  without  completely  restructuring  the  A-matrix.  To  do  so  we 
Introduce  a  vector  0  of  "free”  parameters  which  are  related  to  the 
original  parameter  vector  via  a  matrix  T,  that  is: 

Ap  -  T  0  (11) 

For  the  general  case,  eq.  (7)  and  eq.  (11)  are  to  be  combined: 

AG  -  A  T  A0 

Details  on  estimating  resolution,  standard  errors,  etc.,  for  the  ori¬ 
ginal  variables  using  the  eq.  (11)  transformation  can  be  found  in 
Husebye  and  Hovland  (1982). 

Gravity  inversion  results 

All  the  8  residual  gravity  anomaly  profiles  (Ill  to  X)  shown  in 
Fig.  VII. 10.2  have  been  analysed  using  the  inversion  scheme  outlined 
in  the  previous  section.  For  illustrative  purposes  we  only  display 
results  for  3  of  the  profiles  (Figs.  VII. 10. 3-5),  namely.  III,  VI  and 
IX,  as  these  profiles  are  representative  of  the  northern,  central  and 
southern  segments  of  the  Oslo  Graben.  As  regards  practical  details  of 


the  inversion  analysis  like  initial  and  final  numbers  of  prisms  used, 
number  of  eigenvalues  ignored  in  thesolution,  etc.,  this  information 
is  included  in  the  respective  figure  captions.  Furthermore,  similar 
results  and  details  for  the  other  profiles  can  be  found  in  Wessel 
(1984). 

The  two  striking  features  characteristic  of  all  profiles  are:  i)  the 
small  density  contrast  resolved  and  ii)  the  depth  to  the  causative 
body.  No  higher  value  than  0.06  cm  g-3  for  the  anomalous  body  wave 
were  encountered  during  the  analysis,  while  the  corresponding  excess 
masses  are  concentrated  in  the  depth  interval  15  to  20  km.  The  only 
exception  here  is  profile  111  where  significant  excess  masses  are 
found  down  to  25  km.  The  size  of  the  anomalous  body  gradually 
increases  southward  implying  a  widening  of  the  graben  in  that  direc¬ 
tion  as  Indicated  by  the  surface  geology  as  well.  Furthermore,  the 
very  small  discrepancies  between  observed  and  calculated  gravity  are 
reflected  in  low  relative  RMS-values  of  the  order  of  0.05  mgal. 

Very  recently  Kibsgaard  (1985)  has  studied  fault  mechanisms  of  small 
earthquakes  in  the  northern  part  of  the  Oslo  Graben  (close  to  profile 
III).  The  very  interesting  feature  here  is  that  foci  positions  appear 
to  be  located  on  the  flanks  of  the  anomalous  body  derived  from  the 
gravity  study  (see  Fig.  VII. 10. 6).  A.  tenative,  physical  explanation 
here  is  as  follows;  the  primary  lithosphere  stress  field  in  Fennoscan- 
dia  is  caused  by  ridge  push  forces  originating  in  the  Norwegian  Sea 
(e.g.,  see  Husebye  et  al,  1978)  whose  orientation  in  the  graben  area 
coincide  with  the  compression  axis  reported  by  Kibsgaard  in  his  study. 
The  Interaction  of  this  field  with  loading  stresses  caused  by  the 
"heavy"  intrusive  body  in  the  graben  give  rise  to  stress  amplifica¬ 
tions  being  largest  on  the  lower  periphery  of  the  body  and  ultimately 
becoming  causative  of  the  majority  of  local  earthquakes  observed. 


Extensive  details  on  the  gravity  inversion  investigation  reported 
here,  including  hypotheses  on  Oslo  Graben  taphrogenesis ,  can  be  found 
in  Wessel  (1984)  and  Weasel  and  Husebye  (1985). 

Concluding  remarks 

The  crustal  structure  model  presented  here  and  derived  primarily  on 
the  basis  of  gravity  observations  is  to  our  knowledge  the  most 
detailed  ever  presented  for  the  NORSAR  array  siting  area.  Important, 
work  is  already  under  way  to  compute  the  "seismic  responses"  of  these 
structures  on  crustal  travelling  waves.  Part  of  the  motivations  for 
undertaking  such  work  is  that  preliminary  3-comp,  seismogram  analysis 
results  clearly  demonstrate  a  consistent  migration  off  azimuth  of 
later  arriving  phases. 

There  is  also  another  aspect  of  the  Oslo  Graben  results,  namely,  the 
underlying  dynamic  mechanisms  causing  the  formation  of  the  relatively 
dense  body  in  the  central  crust.  To  improve  our  ability  to  handle  this 
kind  of  problems,  a  cooperative  research  venture  with  Professor  H.-J. 
Neugebauer,  Clausthal,  FRG,  is  in  the  planning  stages.  Tangible  $ 

results  here  so  far  are  a  simplified  modelling  experiment  with  magma 
intrusion  in  the  crust  as  demonstrated  in  Fig.  VII. 10. 7  and  conducted 
by  Dr.  Neugebauer. 


E.S.  Husebye 
P.  Wessel,  now  at  LDGO, 
Columbia  University 
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Fig.  VII. 10.1  The  geology  of  the  Oslo  Graben  and  adjacent  areas. 
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Fig.  VII. 10. 3  Inversion  results  In  terms  of  2i~D  density  contrast  dis¬ 
tribution  obtained  for  profile  III.  The  initial  model  con¬ 
sisted  of  44  prisms,  finally  reduced  to  35  by  combining 
those  having  equal  density  contrast.  Density  contrasts  for 
the  non-hatched  prisms  amounted  at  most  to  ±  0.01  gcm~3  and 
this  was  not  considered  significant.  Observed  (vertical 
bars)  and  calculated  (solid  lines)  gravity  (RMS  *  0.045) 
together  with  the  Moho  contours  are  also  shown. 


Fig.  VII. 10.4  Same  as  Fig.  VII. 10. 3,  but  with  initial  model  consisting  of 
48  prisms,  reduced  to  35.  RM  value  =  0.050. 


Fig.  VII. 10.5  Same  as  Fig.  VII. 10.4,  but  with  Initial  model  consisting  of 
46  prisms,  reduced  to  33.  RMS  value  -  0.049. 


Fig.  VII. 10.6  Projection  of  hypocenters  into  the  profile  IV  cross-section 
It  is  rather  remarkable  that  all  events  are  located  at 
(wrapped  around)  the  lower  part  of  the  gravity  derived  ano¬ 
malous  body.  The  poorest  fit  is  seen  for  event  7  whose  pro¬ 
jection  into  profile  III  gives  a  much  better  fit  vs  the 
location  of  that  anomalous  body. 
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Fig.  VII. 10.7  Synthetic  modelling  of  crustal  intrusions  for  providing  a 
deeper  insight  on  the  mechanism  by  which  the  relatively 
dense  body  in  the  Oslo  Graben  was  formed.  Thick  lines:  boun¬ 
daries  of  intrusions;  thin  lines:  temperature  intervals  of 
200°  (zero  at  surface  and  1200°  at  bottom  or  Moho).  Origi¬ 
nally  the  model  had  one  homogeneous  intrusion,  vertical 
viscosity  contrast  being  5  orders  of  magnitude.  In  the 
beginning  the  intrusion  would  represent  a  density  inversion. 
When  cooling  and  constracting  the  density  contrast  would 
be  positive,  but  located  in  the  lower  crust,  contrasting  our 
gravity  results  for  the  Oslo  Graben.  All  calculations  here 
have  been  made  by  Professor  H.-J.  Neugebauer,  Clausthal, 

FRG. 


